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I— I Abstract 

^ We present a microscopic formalism that extends the traditional formulation of Williams, 

Ericson and Bloch and permits to obtain the transition strengths (TS) of pre- equilibrium nuclear 
reactions directly from their quantum microscopic description. We calculate the TS without 
resorting to the Laplace transform approach and the use of the saddle point approximation. We 
\^ also analyze some problems that may appear in connection with these mathematical tools and 

the Darwin-Fowler approach in this case. 
T-H We show that, analogously to the nuclear densities, the strengths for transitions that change the 

^ exciton number by two or leave it unchanged can be estimated microscopically as convolutions 

of the functions of simpler states. When using the HO basis for the Model Space we obtained 
important departure from the results of the exciton model (EXM), which can partially invalidate 
our previous analysis on the attainment of equilibrium during the PE stage. On the other hand, 
by using constant grid of energies for the sp-basis we were able to reproduce the results of EXM 
quite well in a large range of excitation energies. 

A new model code, TRANSNU, was developed that can be ported to traditional semi- classical 
codes like TNC for nuclear data evaluation. 



1. Introduction 



A common description of the preequilibrium stage of nuclear reactions (PE) uses the exciton 
model (EXM),[T1 [21 [3l H] to analyze nuclear states and particle emissions before the formation of 
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the compound nucleus. In this model the nuclear states are analyzed in terms of their "complexity" 
defined by the number of excited single particle states (sp-states) in comparison with the funda- 
mental state, where the total nuclear excitation is zero and all component nucleons occupy the 
lowest possible energies levels up to a maximum called Fermi level, Ep. 

To describe the excited nuclear states, one initially considers the fundamental state and the 
sp-states that can be excited on it, and call them "holes" when created below Ep and "particles" 
when created above it. In addition, all initial sp-states above Ep are also defined as holes. If at 
a given moment of the evolution of the nuclear system there are "/i" holes "p" particles the total 
number n=p+h of excited sp-states is called the exciton number of the system. For nuclear states 
with excitation energy greater that zero, one assumes that new particles can be created only in 
states previously occupied by holes and, vice-versa, holes can only be created in states previously 
occupied by particles. 

When using the Model Space description though, these simple definitions are expanded to 
consider particles and holes as independent fermion fields, i. e., which can be created or destroyed 
independently. This is achieved by defining a sequence of independent possible levels for each field 
and Ep as an independent phenomenological parameter of the model. Another central aspect is the 
basis adopted for the single particle wave- functions, which in this work is the Hamonic Oscillator 
(H.O.) one.[5j 

The connection with actual physical systems is realized by assuming that the sp-energies 
cannot be greater than a given phenomenological maximum and that particles and holes can only 
be created or destroyed simultaneously, i. e., in pairs particle- hole sp-states (ph-pairs). Then the 
total number of excitons can only vary in steps of ±2. 

From the initial excited state the system is supposed to evolve by increasingly sharing the 
total excitation among the largest possible number of sp-states. The sharing occurs at the mi- 
croscopic level as a consequence of the interaction between sp-states, yielding on the average a 
systematically growing number of excitons up to a maximum, which is the "most probable exciton 
number" at the PE stage. 

On the other hand, the interaction of the excited sp-states is usually defined by a phenomeno- 
logical "residual term" added to the nuclear Hamiltonian and the basis of sp-states can be defined 
self-consistently or phenomenologically. 

The increase of n is counterbalanced by possible PE emission of particles or annihilation 
of a ph-pairs, and after the system reaches some level of "complexity" it is supposed to evolve 
preferrably towards full equilibrium (compound state) instead of emitting more particles at the PE 
stage. 
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In the EXM the rates of transition between nuclear states of increasing complexity are usually 
given in an essentially phenomenological way. In particular, the transition matrix elements can be 
considered as phenomenological constants or simple functions of the nuclear excitation energy, U, 
as exemplified in Ref.[6]. 

In addition, the EXM relies on approximations of statistical nature that may not be very 
precisely defined at the microscopic level, which often makes it difficult to evaluate the importance 
of the details of the microscopic interaction in the description of the PE process. For example, 
the semi-classical formulations^ give the density of states for a given U as a convolution of the 
densities for p and 

rU 

iOph{U)= / ujp{E)LOh{U-E)dE , (1.1) 

JO 

which results from the continuum approximation {CAP) and does not come directly from the 
quantum microscopic description. 

The traditional approach to the Shell Model defines the moments of the Hamiltonian in terms 



of Laplace transforms and their inverse to obtain expressions like (1.1) for the nuclear density and 
in the microscopic description of Ref. [8] a similar approach is followed to obtain the expressions for 
the transition strengths. 

In this work we follow closely the approach of Ref. [8] and present a detailed direct microscopic 
formalism in which one is able to produce the results of the Shell Model without effectively having 
to resort to the Laplace transform and, therefore, without having to rely on the validity of CAP 
and other common approximations of the semi-classical models. 

In this respect the direct microscopic formalism is a more natural and intuitive description 
and it reveals itself as a more appropriate tool to describe, for example, the nuclear transition 
strengths, bringing a complementary view to the usual statistical approach to describe PE states. 

One of the central aspects of this formalism is the proportionality of the degeneracy of a 
given nuclear state (with given total excitation energy, [/, angular momentum, M, and number of 
excitons, p and h), d{U,M,p,h), to the corresponding density of nuclear states 

d{U,M,p,h) (xuj{U,M,p,h), (1.2) 



which permits to connect to the traditional approach, using CAP, and to reinterpret the moments 
of the nuclear Hamiltonian in this limit in terms of convolutions over the excitation energy of the 
moments of less excited states in accordance with the general proposal of Ref. [8]. 



4 



In Sec. 2 we analyze the definition of the nuclear density in the direct microscopic formalism 
in connection with the degeneracy of nuclear levels, the basic definitions of the Shell Model and 
the Darwin-Fowler statistics. 



In Sec. 3 we present the formal definitions of the direct microscopic algebra and present the 



specific ones that apply to the computation of the momenta of the nuclear Hamiltonian. In Sec. 4 



we use the limit of CAP and the usual description of the grand canonical ensemble of Statistical 
Mechanics to make the connection with the Laplace transform, in agreement with the traditional 
approaches. 

Sections [5] and [6] describe the specific application of the direct approach to the evaluation of 
the momenta of the one-body interaction and the residual term of the PE Hamiltonian, respectively, 



and in Sec. 7 we present our results and conclusions. 



2. The nuclear level density 



In the usual "statistical description" of the nucleus inspired by the Shell Model, [3] the grand 
canonical ensemble can be defined by the following relation for the nuclear mass number, 

A = Y,ni, (2.1) 

(i) 

where n, G {0, 1}, are the occupation numbers of the single particle (sp) states associated with the 
corresponding set of sp-levels with energies 

ei = Vie . (2.2) 



The f j are integers and e is a fixed real number, defining the spacing between any two consecutive 
sp-levels of the approximate "equidistant spacing model" . e can also be considered as an average 
spacing of more realistic bases for the sp-states as, e. g., the H.O. basis. The total nuclear energy 
is then given by 

E = Me = Y^UiUie . (2.3) 

If a given physical quantity has discrete eigenvalues defined as a function of a set of integers, 
we will call the latter the set of characteristic integers associated with that quantity. Then, in the 
above case, the set {ui} are the characteristic integers associated with the energies of the sp-levels, 
and the integers {rii} indicate the occupations of the sp-states as either filled with one fermion. 
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= 1, or empty, rii = 0. This specific definition can be useful when dealing with the microscopic 
definition of nuclear states explicitly in terms of the component sp-states. In this case one may 
have the formation of quasi-continuum sequences of nuclear levels, which nonetheless will still be 
univocally related to a set of integers because the bound component sp-states have discrete energies. 

If more than one element of the set {fj, i = 1, .., oo} are equal, the corresponding elements, 

t'ii = • • • = i^ifc e = 1, ••,00} , (2.4) 



define the set of the degenerate sp-states of the energy level 

Ci = i^iie = ■ ■ ■ = i^ike. (2.5) 



Similarly, the number of elements of the set of all configurations of sp-states (nuclear states) asso- 
ciated with a given nuclear energy E is, by definition, the degeneracy of the corresponding nuclear 



level. For given A and E, the nuclear degeneracy is equal to the number of solutions of the Eqs.(2.1 ) 



and (2.3). Therefore, if one considers the microscopic distribution of nuclear states as a function 
of the nuclear energy, the cummulative number of states increase in steps equal to the degeneracy 
of each nuclear level. 

In analogy to sp-states, the degeneracy of a given nuclear level E can be considered as a 
possible characteristic integer associated with E, as the occupation number or "population" of 
configurations associted with each level. 

This number is usually large because there are many possible sets with A non-null elements, 
{riji, • • - jriiA}, satisfying these equations, with values of Ui varying from to A/". The degeneracy 
as a function of E depends on the reference energy used to define the levels of particles and holes. 
For exemple, when the description of sp-states is made in terms of excitons, the particle and hole 
states may have energies conveniently defined with respect to ep and the nuclear energy will be 
equal to the excitation energy. In this case, the ground state has zero energy and all characteristic 
integers associated with the energies of the sp-levels, {fj}, are zero. 

Notice that, in the microscopic description the degeneracy of the sp-levels is defined by the 
specific features of the system hamiltonian, more precisely the symmetries of the nuclear system, 
while the nuclear degeneracy has a more "combinatorial" meaning, in terms of the distribution of 
nucleons into a "pre-defined" structure of sp-states. 

In this context, the nuclear level density at the energy E, p{E), can be intuitively defined 
as the ratio between the nuclear degeneracy and the sp-level spacing e. This is reasonable because 
when the nuclear energy varies from a given value E to its next Enext, corresponding to the variation 
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from the corresponding characteristic integer to its next, the cumulative number of nuclear states 
varies by an amount equal to the degeneracy of the nuclear level at E^ext and, by hypothesis, the 
corresponding variation in the nuclear energy, Enext — E, should be proportional to e, 

Enext - E (X€ . (2.6) 



Now, if one considers the continuum approximation (CAP) and assumes that the nuclear 
energies vary very little from E to Enext , then the assumption of these variations to be proportional 
to e forces it logically to be even smaller. Then in the limit of a set of nuclear states with "very 



large" density Eq.(2.6) becomes, 

Enext -E ^ dE ^ e (2.7) 

and the corresponding variation of the cumulative number of states would be a function of the 
energy given by 

dM ^Mnext-M = Vk{E), (2.8) 

where T>k{E) is the nuclear degeneracy at the excitation E. 

Then, the nuclear density p{E) becomes approximately. 

The same result can be obtained using a more formal argument, by considering e constant in 



(2.3) and applying the corresponding definition of the statistical Shell Model formalism using the 
Darwin- Fowler method. [3] The generating function for the corresponding grand canonical ensemble 
is given by the expression 

f{x,y) = J\{l + xy''^) = J\{l + Xi) , (2.10) 

i i 

where the last simple form takes into account the fact that x is indiferentiated for the various 
component sp-states as it accounts for the number of "relevant" sets of sp-states independently 
of their specific characteristics, i. e., x is the same for all sp-states. Both x and y have physical 
meaning under the statistical description of the grand canonical ensemble, for many-body system 
with variable number of "bodies" and variable energy, but x has a more strictly combinatorial 
meaning while y is related with the probability distribution associated with the various microscopic 
systems of the ensemble. 
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Then, the nuclear level density can be directly defined as an adequate pole of the generating 
function divided by e,[3] 



while the generating function can be rewritten as 



f{x,y)dxdy 



(2.11) 



f{x, y) = l + x^y^J +x'^ y\-ji^-j2) _^ 



^ ^ y^'^J^^ ti'ja; _|_ . . . ^ 



(2.12) 



A y{^j\^ Vvja) 

(i) (jij2) 

which, therefore, describes all nuclear systems with all possible "mass numbers" (relevant sp-states 
in each microstate) and energies (total energy of the microstate). 

In other words, for nuclear systems each configuration of sp-states is also a microstate of the 
canonical ensemble with fixed mass number and temperature, [9] and the term proportional to x^ 
is the sum over all possible configurations with fixed nuclear mass A and variable energy. 

For each nuclear level, E'fc, corresponds usually many different configurations of sp-states and 
to each nuclear mass A a term, y^-, is defined in Eg. ( 2.12"| ) as follows. 



(2.13) 



which can be rewritten in terms of the degeneracies for the various nuclear levels, as follows 

(2.14) 



where 



(2.15) 



i=l 



and Dk = T^k{Ek, A) is the degeneracy of the nuclear level E^, for a given nuclear mass number A. 
Taking Eg. ( 2.14| ) into Eg. (2.12) and using the definition (2.11) for the nuclear level density 

yields 



p{Ek,A) = -VkiEk,A) 



(2.16) 



where usually the dependence in A is not explicitly written. Therefore, the definition in Eg. (2. 11) 
is coherent with the "intuitive" description given at the beginning of this section. 
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Equation (2.12) can be rewritten as 

fix, y) = l + xYPki{Eku + 
(fei) 

x^YPk2{Eu2, 2)?/-^-(2) + • • • + + • • •+ (2.17) 

m (k) 

where different indices liave been used for each term to reinforce the fact that the corresponding 
nuclear levels may not be the same. In these expressions the sum over k is equivalent to the sum 
over Ek, then f{x,y) can be rewritten as 

f{x,y)= ^P(AC/)xV/^ (2.18) 

{A,U) 

and also as a sum over individual configurations, with all degneracies are equal one, 

f{x,y)= J^xV/^. (2.19) 

(con/) 



3. Formal definitions of the direct microscopic approach 

To describe a given nuclear excited state, with variable mass number A and energy E, one 
may assume that it belongs to the corresponding grand canonical distribution defined over the 



nuclear configurations obtained as the solutions of the Eqs.(2.1) and (2.3) 



In this case one may always assume that the elements of the set of characteristic integers of 



the energies of sp-states, {ui} in (2.3), are ordered according to increasing values as a function of i 
and, by definition, there are only A nonzero elements in the set {rii}. 

For a given nuclear mass A the nuclear state, , can be represented by a product of sp- 
states, which are usually approximated by self- consistent quasi-particle states, and can be written 
as linear combinations of the corresponding complete set of eigenvectors of the sp-hamiltonian, 

A A 
i=l *=l(ki) 

= E- • -E^i^i- • • . (3.1) 

(ki) (kA) 

which can be rewritten as 

I^a) = j;C,e^^^*/^|W,) (3.2) 
(j) 
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where j is a characteristic integer for the nuclear energy, corresponding to a specific sequence of 
the enumerable set of components of In this case the notation (j) is used to designate the 
corresponding degenerate set of nuclear configurations with energy Ej . 

The grand canonical ensemble vector corresponding to the various nuclear states \iPa) can 



also be written[TO] directly in terms of the occupation numbers of all sp-states defined in (2.3) as 

(ni,---,"oo) 

and if one defines the symbol "(ki- • -k^)" for the set of all configurations containing A and only A 
sp-states with non null occupations as 



(kr • -k^) = {{ki, • • ■,]^a}/- n-^i = ■ ■ ■ = n^A = l}' foi^ given A, 



(3.4) 



where ki, • • •, k^ are also supposed to be ordered by increasing values, then ( |3.3[ ) can be written as 
IV') = X] |m,---,"fci,---,f^fcA,---,"-oo) = ^ ^ |(kr • -k^)) = ^IV'a) (3.5) 

(•■■,ki,-,kA,-) (^)(ki,-,kA) {A) 

The Fock space operator of the excitons, "particles" (p) and "holes" (/i), of an arbitrary 
nuclear system, with arbitrary A, that describes the corresponding grand canonical ensemble of 



nuclear states, with expected values on the states |^) given by (2.10), can be defined in the case 
of two fermion fields, "p" and "/i", as[S] 



111/ iiu 

If one considers initially the simpler one fermion expression 

-^0 = + Xpia\ai 



(3.6) 



(3.7) 



then it can be rewritten as 



oo 

n[ 

1=0 



aia'i 



(sl<s2) 



n 

l'{ll.l2.l3)^ 

LV 7^(s1,s2) > 



+ E 



n 



mail 



ana'ii 



-A 



.t 



••I 0/20^2 I- •• 



(3. 



XpslXps2 (3.9) 
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+ X] (• • ■)XpslXps2Xps2. H (3.10) 

(sl<s2<s3) 



oo IN 



{Af=0){s^r) Vj=l 



where, for fixed and configuration (sat), n(s7v) is given by 
n(s7v) = n ■ ■ ■ • • • (oiivaw) • • • ("Ijv"sA^) • • • ("«(iV+i)4(7V+i)) ■ ■ ■ (3-12) 

/(ll.---,l]v,liV+l)N 
V 7^(sl,.--,s]V) / 

where there are terms of the type ^ajj a^j ^ and infinite terms of the type ^o/ j a|j ^ . 
Then it is clear that 

n(s^)|(kr--kA)) =5^,A5(sr--sjv|kr--kA)|(kr--kA)) (3.13) 

and therefore 

Fo\ij) = FoY,\^a) = (3.14) 

{A) 

EEfn^P% )n(s7v)EEl(ki---l^^)) = (3.15) 

{N){st,)\i=l J (A)(k^) 

E E ( n^p^j'^^'^'^(^i---^^i^i---^^)i(^i---^^)) = (3-16) 



E E (n^p«. )i(^i---^^))=Ek^) (3-17) 

W(^i-siv) Vj = l / W 



which, by comparison with (3.5) shows that Fq projects the grand canonical ensemble vector 



of Eq.(3.3) into another linear combination of its various components, where the coefficients have 



changed from 1 to 

llxps}j , (3.18) 

which are proportional (for given N) to the canonical ensemble probabilities of the component 
configuration |si- • -sat) .[9j 
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Eq. (3.17) also implies that for two given independent grand canonical vectors 



|1) =EEl(ki---k^i)) a^^d |2) = ^ j;|(lr--l^,)) (3.19) 



results, 



N 



(1 1 Fol 2) = ( 1| J] J] Hxps, ]\{sv -sn)) (3.20) 
W(^iv) \j=i 



= EE ( n^p^. ) HsN\msN\2) (3.21) 

Let O be an operator on the Fock space of the sp-states, which can modify the nuclear 
configuration. Then one can write the corresponding transition strength, defined as the square of 
the transition moment (1|0|2) summed over all possible transitions, as|llj 

So= E |OP = El^l|0|2)l' = E^l|0|2)(2|0^|l> 

(trans.) (12) (12) 

due to possible variation in the number of excitons and nuclear excitation Sq must be redefined in 
the framework of the grand canonical ensemble using the grand canonical distribution, which can 
be written schematically as 



|0|2 = (FoOF^Ot) = ^(1 1 FoOF^Ot 1 1) = ^(1 1 Fo| 2) (2 1 0\ 3) (3 1 F^| 4) (4 1 Ot 1 1) (3.22) 
= EE ( n^^'^: ) i^N\0\rM)Y.Yl friv.) {rM\0^\sN) (3.23) 



In the case of two Fermion fields (e. g., "particles" and "holes") Fq is given by Eq. (3.6) 



Fo = FpFh = n + xp^.alaf,j [| [bf^bl + Xhf,b%j , (3.24) 

fj, V 

which can be rewritten as 

oo / p \ 

^v = Y^\Vi.^vsM{^v) (3-25) 

and 

oo I ^ \ 

i^/^ = EE n^'^-: (3.26) 
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and, as we saw in (3.13), due to the properties of the single particle fermion operators one can also 
identify n(sp) and Il{rh) with the components of the respective projection operator, 



n(sp) = \splsp\ and U{rh) = \rhlrh\ 



(3.27) 



and 



U{sp)U{rh) = \sprhlsprh\. 



(3.28) 



Then the analogous of (3.20) becomes 



{l\FpFh\2) = {l\^Yl in ^Ps,^hrAs{sprf,\2)\{sprh)) 

{ph){sprh) \k,j=l ) 



(3.29) 



ph 



2Z 2Z n ^v^k^hr^ \b{sprh\\)b[sprh\1) 

{ph){sprh) \k,i=l I 



and 



p'h' 



(11^X12) = (1| \\{VhrA5{sp'rh'\2)\{sp, 



rh' 



{p'h'){s,r^,) \l,i=l 



(3.30) 



(3.31) 



p'h' 



{ Yl^'psi^'hrA^{sp'rh'\'i^)S{sp'rh'\2) 



{p'h'){s,r^,) \l,i=l 



and the analogous of (3.22) and (3.23) are 



FpFf,nF;F'^0A = Y,{MFpFh^K^LO^\'^) 



(1) 



(3.32) 



(3.33) 



^(l|F,F,|2)(2|J]|3)(3|F;F;:|4)(4|Ot|l) 
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(3.34) 



ph \ / p'h' 

Yl^pskXhrA{sprh\^\sp'rh>)J2 Yj \ Yl^'psi^'hn\{sp'rh'\0''\sprh) (3.35) 



{ph){sprh) \k,i=l 



{p'h'){3^,ry) \l,i=l 
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which can be rewritten as 

(ph p'h' \ 
n n ^P^kXhvjX'p.^x'h^^ {spTh I n| SpTh'Xsp'rh' I O^Sprh) (3.36) 
k,j = ll,i=l I 

Now we make the usual change of variables that defines the explicit connection with the 
microscopic statistical parameters of the sp-states and also brings the possibility of the Laplace 
transform interpretation 



( — Be,-. — 7nir; ) 



xL=y'e(-^'^'-»-V-'^^) (3.37) 



then (|3.36|) can be rewritten as 

ph p'h' 



Kp'h'J '^V'V^ 

{spTh I 0.\ Spirh'\spirh' 1 0^\sprh) (3.38) 



or 



= J^xV'^'yS"^' ^-'^''-'''-^''''-''''\sprn\n\sp,rt,,lsp.rn>\0^\sprn) (3.39) 

\p'h'l \s^,r^,J 

where 

ph ph 

U = Y esk+ and M = ^ nis, + m^^ (3.40) 
(fc,j=i) (fc>j=i) 
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p'h' p'h' 

U' = ^ esi + and M' = ^ + nir, (3.41) 

{l,i=i) {l,i=i) 

As a shorthand practical notation that includes the essential features of the above expressions 
one may define 

(FpFh^F'j^FiO^) = {VL\0^) = {sprhMsp>rh,lsp,rh'\0^\sprh) (3.42) 

(12)(!7M) 



where 



^^e[^^](spr^|0| vr^'Xv^/^' W (3.43) 
(1) (2) 



E = E ^"^'^V'y"' (3.44) 



represents the sum over all possible numbers of excitons and 



(3.46) 



(2) ( 'v-h ) 

represents the sum over all configurations for a given exciton number and 

QpM] ^ ^-pU-iM-l3'U'-i'M' 

which is the non normalized gand canonical distribution function. 
At last one may just drop the s's and r's and write 

(Ji|ot) = EE^''''''(^'^l^l^''^'X^^'^'|O^I^'^)' (3-47) 

(1) (2) 

which now has a precise meaning, where \p, h) , \p', h') represent the possible configurations for 
given exciton numbers p,h and p' ,h' . 

In the case of the simple expected values of Q, the expressions are totally analogous, 

{FpFhn) = {n)= E {sprhmsprh) = J2J2^^'^''^^'p''h\^\^p''h) (3-48) 

(12)(C7M) (1) (2) 



with the two sums given by 



E = E^'2/' E = E (3-49) 

(1) iph) (2) (spru) 
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and 



(3.50) 



Again one may drop the s's and r's to obtain the simphfied expression 

(1) {ph) 



(3.51) 



Now, it is clear that the appHcation of CAP on Eq. (3.44), for example, will produce a very 
large number of levels per unit energy on ^2 ^^^d permit the approximate replacement of the sum 
by an integral. The details of this procedure, its interpretation and consequences are analyzed in 
the next section. 



4. The connection with the Laplace transform 



From the definitions of Sec. 2 it results that the nuclear excitation, U , is a parameter that 
varies between two finite extremes 



Emin ^ U En 



(4.1) 



where it takes a sequence of discrete values with degeneracy 'D{A,U ,M), as defined by Eqs.(2.1) 



and (2.3), with A interpreted as the total number of excitons, 

A = n = p + h. 



(4.2) 



Similarly the total nuclear momentum also varies in a stepwise manner between M^in and 
Mjnax and any additive quantum number of the total system can be treated similarly. 



As it is shown in Eq.(2.17), in expressions like (3.43) or (3.47) for given numbers {p,h), there 
is in general a subset of the configurations {spTh) for which the quantum numbers {U,M) take the 
same values and, by definition, the number of elements of this subset is equal to the degeneracy of 
the corresponding nuclear state. 



{isp)virh)i} = {{si, 



n, • • • 



rh}i, i £{!,■■ -MA U,M)}} . 



(4.3) 



Then, e. g., one can rewrite Eq.(3.48) as 
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(12) 

Mma^ Emax V{A,U,M) 

(1) M=M™,„ U=Emi-a * = 1 



(4.4) 



From the phenomenological calculations with nuclear level densities [12] one knows that the 
cumulative number of nuclear states can be very high even for not very high excitations. Then, 



it is clear that the use of CAP to replace the sum over U by an integral in (4.4) is a reasonable 



procedure, although approximate, and one can explore this possibility using an ad hoc definition 



of nuclear density, inspired by the analysis of Sec. 2 



Note that for each U in the sum in the RHS of (4.4) there are T>{A,U ,M) states with the same 
energy and angular momentum, for which e^^*^! has the same value. Then, if 5U={U — Uprev) is 
the variation of the nuclear excitation between its present and "previous" value, the corresponding 



approximate nuclear density for each index i on the last sum of the RHS of (4.4) will be 



uj{A, U, M) V{A, U, M)/6U = constant, 



and 



V{A,U,M) 



V{A,U,M)= {l)^uj{A,U,M)6U= uj{A,U, M){l)dU. 

(^i = l"j Uprev 



(4.5) 



Then, replacing "(1)" by an arbitrary integrable function "(...)" and summing over all U gives. 



Ema:. V{A,U,M) 

E E (■■■ 

{U = Ermr,) (i=l) 



U = E„ 



Uj{A,U,M){...)dU. 



(4.6) 



and the sum over configurations with given a number of excitons in (4.4) becomes 



Y Y Y 

M=Mrmn U=Emin « = ! 



-I3U 



u}{A,U,M)e~^^ndU 



(4.7) 



JM=M^i„) 



where the definition (3.50) for e^^^^l was used. 



Now for the nuclear excitation energies the minimum is the ground state corresponding to 
Emin = and the maximum is unbounded and one can take,|T] 



Emax ~ CO) with good approximation, 



(4i 
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then (4.7) becomes the Laplace transform of the part of the integrand inside the parenthesis. 



More generahy, one can write (4.4) as 



{Fj^F^n) = Y. Yl E Y.^^'^'^^daiP, h, U, M){n)^{p, h, U, M) (4.9) 

(1) (M=M^i„)(t/=0) («) 

where a indicates ah configurations for which "(p/i[/M| il|p/iC/M)" has the same value, i. e., 
the configurations degenerated with respect to the action of VL or the observation of the physical 
quantity represented by fi, and ^^da{p,h,U,M)" is the corresponding degeneracy of states, and 

Yda{A,U,M) =V{A,U,M) . (4.10) 

(a) 



Then, the expression corresponding to (4.7) is 

ax /» CO 



{FpFhn)^Yl E / '^f^e-^''-'^J^a;,(p,/i,[7,M)(f7)„(p,/i,[7,M) (4.11) 

(1) (Af=Af„,„)-^° (a) 

where the nuclear density is approximately defined as 

cOaiA, U, M) da{p, h, U, M) /5U = constant, 
for each a, and one can rewrite 

d„(p, h, U, M) K Wa(A, U, M)x{U - Uprev) ~ i^a{A, U, M)6U , (4.12) 
which corresponds, for given number of excitons (p, h), to 

^ / dUe-^''-^''Yu;^ip,h,U,M){n)^ip,h,U,M) . (4.13) 

(M=Mrninr° (a) 

Therefore, by definition, 

e-^^'^Yua{p,h,U,M){p h UM\VL\p h UM)^ \ , (4.14) 

(Af=M„,„) (a) I 



where the symbol C{X} indicates the Laplace transform of {X). 
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In particular, in the case of the unitary operator, Q=l , all expected values are equal "1" 
and {FpFh) becomes essentially the Laplace transform of the state density, then 

e-^'''u{p,h,U,M)^C-\{FpFh)) = C-^\^llil + Xp^)l[il+Xh.)] . (4.15) 



In Appendix I it is shown that this result is equivalent to the traditional one of (2.11 ) because 



the RHS of Eq.(2.11 ) can be interpreted, using CAP, as the inverse Laplace transform of the grand 
canonical generating function. [H |H] 

Therefore, the microscopic formalism yields the expression of the inverse Laplace transform 
of the expected values of the interacting operators without having to actually evaluate it. 



4.1 Convolutions and discrete functions of the configurations 



Regarding the idea of nuclear degeneracy, one notices that in (4.4) all the expected val- 
ues "((sp)j(r/i) ■ I r^l (sp) -(r/i)^)" correspond to configurations with energy U and total degeneracy 
P(yl,C/,M), which we will designate by {u) 

{u) = {{spurri);i = l,V{A,U,M)] , (4.16) 

with A=n=p+h. 

If the operator Q. describes the measurement of the spin of the nuclear state, for example, 
it will have in general a sequence of different discrete values for the different elements of {u) and 
the expected values, (il), will also be degenerated, i. e., in general there will be more than one 
configuration for each value of the total spin. Then, one may write 

V{A,U,M) V{A,U,M) Pma. 

M (i=l) (i=l) /9=i 

{conf) (spins) 

where {/3 = l,/3max} is a sequence of integers in biuni vocal (one to one) correspondence with the 
discrete set of "spin values" |(r2)^|^ and di3{A,U,M is the "spin degeneracy" satisfying, 

YdM^U,M)=V{A,U,M) (4.18) 
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This idea can be straightforwardly generahzed to the "measurement" of a quantity that 
changes the total number of excitons, A=n={p-\-h), for example: Q=Y^^^-^Oaaaaa, where "a" 



indicates all possible exciton states of "particle" type. Then, the expected values in (4.4) would 



select all configurations in (u), Eq. (4.16), that have one sp-state "a" in them. 

As we have seen, the number of such configurations define their "degeneracy", for which 
the destruction of an sp-state "a" has non null expected value. If this number is designated 
by da{A,U ,M), then the resulting configuration after the "observation" of Oq,, the intermediary 
state, would have the same degeneracy, but the corresponding grid of nuclear energies and angular 
momenta would be displaced by and ma respctively, 

daiA -1,U -ea,M - ma) = daiA, U, M) . (4.19) 

In this context the introduction of CAP naturally brings the idea of convolution between 
states of different levels of complexity as it transforms the sums of discrete sets of values times the 



respective degeneracies, Eq. (4.17), into integrals involving the corresponding nuclear densities. 



The sums and integrals have their ranges defined by the idea of available states, i. e. the 



states that give non zero expected value for the measured quantity. For example, in (4.13) the 
density uJa{A,U ,M) corresponds to available states associated with the operation or measurement 
of 0, in the subset of the grand canonical ensemble with given A, U and M. 

On the other hand, the set of configurations for each a is not directly connected with the 
nuclear excitation energy, but if the density of states is high all configurations with energies between 
two given energies, U and Uprev, can be considered as having energy U and the variation over a 
can be associated with densities corresponding to U. 

In particular, when one sums over an sp-state "a" the previous discussion shows that it 
could have some degeneracy and, therefore, it is not equivalent to the sum over 6^. This idea is 



further developed in Appendix II where the sum over is then reduced to a convolution of nuclear 
densities, i. e., the sp-energies have an associated nuclear density that, in the CAP limit, naturally 
defines a convolution with the nuclear density present in the expected values. 

Notice that the sum over the total nuclear spin projection M, and the corresponding sums 
over single particle angular momenta, ma, cannot be transformed in the same way as U and e^, 
using CAP, because the interval between sucessive values of M or is never smaller than 0.5 and 
therefore it cannot be considered as an infinitesimal even for a large number of configurations. 
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Next we present a brief application of the above formalism for the calculation of the momenta 
of the operator of one body transitions to illustrate the main techniques to evaluate the transition 
strengths (TS) of the PE nuclear Hamiltonian. 

5. Brief discussion about the one body operator 



This section presents a detailed calculation of the transition strengths of a generic one body 
operator using the direct microscopic formalism. The results are given either as discrete sums of 
the microscopic approach or, using CAP, as approximate convolutions involving the nuclear density 
of the intermediary state and the densities of sp-states. 

The one body operator 

O = ^ Oa(3aiai3 , (5.1) 

(a G p\ 

\0epJ 

could describe a physical process involving a single particle transition like a 7 transition or the 
particle emission into the continuum in a direct statistical process. [8j 

When the "particles" are redefined as excitons and the Model Space also includes "holes" it 
becomes necessary to take the contribution of the latter into account. In the EXM the particles 
and holes of the Model Space are considered as independent and complementary fields, in which the 
creation of one exciton can also be described as the anihilation of the complementary one. Then, 
the description of this pair of independent and complementary fields is achieved by replacing 

ap by {ap + b^^) and by (a^ + ha) , (5.2) 



then (5.1) becomes. 



O = ^ Oai3aiai3 + ^ Oai3aib^i3 + ^ Oapbaap + ^ Oapbabl . (5.3) 

(aep\ (aep\ faeh\ fa£h\ 

\pep) \peh) \p&p) Kp&h) 

This expression is a little more complicated than it looks and some epistemological clarifica- 



tion is presented in Appendix III, where it is shown that (5.3) can be rewritten as 



O = ^ Oo,pa)^ai3 + ^ Oaolaa + ^ Oafsa'lb^p + ^ Oapbadp 

( a£-p \ (oefilled-p) (ae-p\ (aetillBd-h\ 

V/3 g f illed-p^ ^ \/3e-hJ I /3 £ f illed-p^ 

(a) (b) (c) (d) 



bibo^Oo^a + Oo (5.4) 

("^~^) (oefilled-h.) 
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(e) 



(f) 



(g) 



where, 



Oo= ^ Oaa= ^ (a|OQ/3(5a,/3|/3) 
(aeall-h) (a,/3eall-/i) 



Then, if one considers, for example, the number of particles operator, N, in (5.1), all Oaf3 



will be one and Eq. (5.4) gives 



{N)= E b^-ba+ E (1). 

(aGfilled-p) (q e f illed-/i) (aeall-h) 



(5.5) 



Then, using the letter g for the maximum number of "p" excitons and b for the maximum number 
of "/i" excitons in the Model Space, results 



{ph\N\ph) =p - h + b 



(5.6) 



which is equal to the total number of already excited "p" sp-states ( "filled" p states) plus the total 
number of "/i" sp-states that have not yet been excited ("empty" h states). 



In Appendix III it is shown that although the above result is true in the strictly combinatorial 



point of view, the correct physical interpretation of (5.5) is {N)=g, which is in agreement with the 



original idea expressed in (5.1) for Oa/s equal to 5(q,^). 



5.1 Transitions that change the number of excitons 



The terms of (5.4) that change the number of particles or holes of the initial state, \ph) , in 
comparison with the intermediary states, \p'h'), are (c) and (d), corresponding to the description 
of pair creation and destruction in the Model Space. 

The corresponding transition strengths can be calculated using the general algebra developed 

m 



Sec. 3 For example, for the term (c) one has 



c c 



E E I P (P^ I ^fl I P I I P^) 

(12)(l/M)(a/3) 



(5.7) 



and initially one defines all the constraints for (ph) and {p'h') determined by the sp-operators 



{l,2,UM){al3) 



(5. 
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factorize the terms of e^'^'^l selected by the destruction operators, 

= EEio./3i^(^) E ey'Mp-i\mh-i\h') (5.9) 

and use CAP to obtain, 

dd' 

EViy)^(/^-i,/.') E e[-/^^-/^^'-^^^-^^^']5(C/ - 5')5(Af - Ox5:|0.,P , (5.10) 

where 

S = U — ea — ] i = M — niQ, — ; 5' = C/' + ea + and t' = M' + iiIq + . (5.11) 



Now one notices that the sum over the intermediary states reduces to only one term corre- 
sponding to the conditions between configurations 5{p — l\p') and 6{h — l\h'), then 



or more explicitly 



(c|ct) « J]J]e[^*^]^0„;3d(p- l/a,/i- l^/3,5,t) 

(1) {S,t) {al3) 



(c|ct) « ^ J^el^*^] ^0«^d(p -a,h-/3,S, t) 

(1) {S,t) (al3) 



(5.12) 



(5.13) 



In Appendix III it is shown that (5.13) can be rewritten as convolutions between nuclear 



densities and the densities of sp-states 



c c 



CAP 



(1) (S,t) 

where e=(ei+e2), m=(mQ, +111/3) and 



^^g[t/A/] fde^u;{p-l,h-l,U -€,M -m)0,{l,l,e,m) , 

(1) (S,t) (m) 



(5.14) 



0,(1, l,e,m) = j de2^Oai3(^il,0,€ - €2, m - in/3)a;(0, l,e2, 



(5.15) 



and the approximate definition of nuclear density of Sec. 4 was used. 

Similar expressions are valid for the other expected values and transition strengths. For term 



(d) in (5.4) the procedure is analogous, 



(d|(i+) = ^ ^{ph\baafj\ph'1ph'\a^b^\ph) 

(l,2,r/M)(«/3) 



(5.16) 
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E Y.^{pmKh^cc)5{p' = /3)6{h' = a)Sip\p' - l)dih\h' - 1) (5.17) 

(12)([/M)(a/3) 



= EE(^) E ei'%ip\p'-mh\h'-i) (5.18) 

(1) (ajS) \ ^ / (p'=p,h'=c\ 

and using the continuum approximation results, 

dd' 

{d\d+) Y.Y.^^''''%,p'-^Ah,h'-i)m - s')5{M - t')EE(i) ' (5-19) 

(1) {St) {a0){j,k) 



where 



then 



S = U + ea + ei3 and i = M + iiIq, + , (5.20) 

S' = U' -ea- and i' = M' - nia - , (5.21) 



{d\S) « EE^'"'"'^ E^«/3^(^^ + 1= /3, + 1= a, 5, i) , (5.22) 

(1) (S,t) (a/3) 

or more expUcitly 

{d\S) « EE^'^""'^ E^"/^'^^^' + /3, /i + a, [/ + + e;3, M + + m^) , (5.23) 

(1) (5,t) (a^) 

where the sum over (a,/3) can be approximated as 

^Oapd{p + P,h + a, C/(^a^/3) + Ca + M(^c.^^) + ma + my?) (5.24) 

(a/3) 



Ooow(p, /i, M) - jde^uj{p - l,h - 1,U - €,M - in)0,(l, 1, e,m) , (5.25) 



CAP 

(m) 

with 



Ooo = J deid€2 E Oa/3^^(l,0,ei,mQ,)a;(0, l,e2,m^) , (5.26) 



(m„m^) 



and 



0,(l,l,e,m) = j de2^OQ,/3a;(l,0,e-e2,m-ni/3)w(0, l,e2,ni/3) , 



The transition strength for the term (a) in (5.4) is given by 



(a|a+) = ^ '^Oai3{ph\a^ai3\ph'\ph'\a'l^aa\ph) 

(12)(!7M)(a/3) 



^ ^ 0^^6{p=am6{p'=P^a)6{p\p')6{h\h') 

(12)(C/Af){Q/3)(a^/3) 



then 



where 



(1) (a/3) 



^ e[^*l<5(p|p')5(/i|/i') 

/ p-l^a(i,h \ 



CAP 



dd' 



a a ' 



5;^e[^*^]5(py)5(,,,,)5([/ - S')6{M - t')E0a/3E(l) 

(1) (5t) (a/3) (j,fc) 



5 = [7 - + e/3 , t = M - ni„ + m/3 , 5 ' = C/' + e„ - e/3 and t ' = M' + m, 



yielding 



(1) (5,t) (a/3) 



1 + l^a = P,h,S, t) 



or more expUcitly 



(a|at) « ^^e[^^^1^0«/3(i(p-a + /3,/i,;7-e„ + e^,M-m<, + m/3) 

(1) {S,t) (a/3) 



which can be rewritten as (see Appendix II) 



f 

"I") / deiV'Oaw(p- t/-ei,M-ma)w(l,0, ei,mQ,) 
Jo r~ 

(iHq) 

- / de (2, 0, e, m)a;(p - 2, h,U - e, M - m) , 

Jo , N 
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where 

Oa=[ [ de2Y,^{l,0,e2,mfs)Oap | , (5.37) 

e=(ei+e2), m=(mQ +111/3) and 

O,(2,0,e,m) = / (ie2 V]OQ,/3tj(l, 0, e - €2, m - m/3)a;(l, 0, e2,in^) , (5.38) 
Jo I ^ 



and the transition strength for the term (e) in (5.4) is given by 



e e 



-) = {ph\hlho\Tlh;y}^\hthfi\vh) , (5.39) 



or 



Y.^e^''^'^Y.'^a^'^^P,h-l + l^P = a,S,t) , (5.40) 

(1) (5,i) (a/3) 



(e|et) ^J^J^^el^*^] 0«;3(i(iJ,/i + a -/?,[/ + e„-e^,M + m^-m^j) . (5.41) 

(1) {S,t) {a/3) 



which can be rewritten as 



e e' 



/ deiy^Of3d{p,h-l,U-e2,M-mfs)co{0,l,e2,rrif}) (5.42) 
JO 

(ma) 

/ deVOe(0,2,e,ni)(i(p,/i-2,f7-e,M-m) , (5.43) 

Jo , ^ 



where 



0/3= /'(iei^^^(l,0,ei,m„)O«/3 I , (5.44) 
\ K) / 

£=(£1+62), in=(nio+ni_g) and 

Oe(0,2,e,m)= / dei V]OQ,/3a;(0, 1, e - ei,m - ma)w(0, 1, ei,mQ,) , (5.45) 
Jo I . 

(ma) 

Similar results are obtained for the expressions involving the crossed products of the terms 
(6), (/) and (g), For example the product i^\g^) yields, 

(^b^) = E Y.^c.^Ol{vh\aW\^U)<^K\vh> , (5.46) 
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Y,OaaO*o6{p = a)6{p' = a)6{p\p')5{h\h') (5.47) 

(12)(l/A/) (a) 



EEO^M'-^] E ei'^%ip\p')5ih\h') , (5.48) 
xy 



(1) (a) \ " / fp'=^.'''=-) 



then 



dd' 

CAP 



where 



giving 



(1) {St) (a) {j,k) 



S = U , t = M , S' = U' and t=M', (5.50) 



(bW) ~ Y.Y.^^^^'^120'^<^Ood(.P= a,h,Uc.,M^) , (5.51) 

(1) {S,t) (a) 



or more expHcitly 



(bW) ~ EE^^^^^'E^"-'^o4p - a + a, /i, [/ - + e„, M - ni„ + ni„) . (5.52) 

(1) (5,t) (a) 

The total of the nonvanishing terms of the sum over (a) is equal to the number of excited 
"particle" sp-states in \p'h') or \ph) , i. e. "p", then if Oaa is approximately the same for all (a), 



Eq. (5.51) reduces to 



{b\g^) ~ J2J2^^^^^0aaO*opd{p= a, h, Ua, M«) (5.53) 

(1) {s,t) 



6. The pre-equilibrium nuclear Hamiltonian 



This section presents the analysis of the momenta of two-body operators, in which the tran- 
sition strengths may describe the pre-equilibrium nuclear transitions such as those in the multi-step 
formalisms and semi-classical models. [8] 
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To distinguish between the "mean-field" and "residual" components of the total interaction 
one may add a one-body kinetic energy term to the total interacting potential and analyze the 
resulting moments of the Hamiltonian. Those invariant with respect to one body transitions cor- 
respond to the Hartree-Fock mean- field description of the nucleus, without particle emission and 
with fixed nuclear excitation, the remaining ones define the residual interaction responsible for the 
pre-equilibrium transitions. [T3] 

The full Hamiltonian that describes the interacting field of sp-states that compounds the 
nuclear system during the pre-equilibrium stage can be defined as|8] 

H = ^ ta/3aiai3 + ^ ^ Vais-ysata'pasay, (6.1) 



which, after normal ordering the single particle operators for particles and holes, see Eq.(5.3), 
becomes 

H = Eo + Hhf + Vres (6.2) 



where the ground state energy is 



Eq = ^ taa + ^ Vapap-, (6.3) 



2 

(aSp) (a/3ep) 

and the Hartree-Fock Hamiltonian is 

Hhf = ^ Kpa^ap + ^ Kphaap + ^ Kpo^b'^ - ^ Kph^ha- (6.4) 

The condition of diagonalization of H^f (invariable nuclear excitation) is then equivalent to 
the diagonalization of h^p which is expressed by 

hap = ^ap^a-, (6.5) 

where is the single particle energy of the sp-state |q) and the terms involving 60,0/3 and a^b~p have 
vanished because particles and holes are supposed to be independent fields and, therefore, cannot 



be associated with the same sp-state (see Appendix HI for detailed comparison of the Model Space 



and physical interpretations). Equation (6.5) corresponds to the idea that in the Hartree-Fock 
description of many-body systems composed by "particles" and "holes" there is no creation or 
destruction of particle- hole pairs (ph-pairs), but only propagation of sp-states. [T3j 
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Therefore, the consideration of possible change in the number of excitons is restricted to the 



action of the residual interaction, Vres in (6.2), which is introduced as a perturbation over the 



mean-field description and it is not necessarily self- consistent as the Hartree-Fock part usually is. 



Using (6.5) Hhf becomes 



Hhf = ^ eaa^Qa - ^ e/3&^^/3 • (6-6) 
(aep) (/3e/i) 

and if the energies of "particles" are defined relative to ei;' as 

ep = €a- ep (6.7) 

and for "holes" 

eh = ep - ep, (6.8) 

then the nuclear excitation energy becomes the simple sum of the p and h energies, 

C/=5^(6p + e,). (6.9) 

Strictly speaking this definition of the nuclear excitation only makes sense if the numbers 
of particles and holes are the same, otherwise it will add an arbitrary phenomenological term 
proportional to ep. This term can be interpreted as the total excitation at the origin of time, but 
it is not important in the present discussion. 



The last term in the hamiltonian (6.2) is given by 

(q/37<5) 



Vres= ^ (^Vai3^sata^bjb:^^ + ^Vai3^sata'^bja^ + ^Vai3-y5atbjb1^bi3 



(a) (b) (c) 

^Vai3jsaaa'^asa^ - Vajs-ysatbjbfsa^ + |Fq,/3^56| 6+ 60,6/3+ 

(d) (e) (f) 

^Vap-rsatb/sasa^ + ^Vais-ysbjbabisa^ + ^Vap-ysbabpaga^^ , (6.10) 

(S) (h) (i) 



where the various letters below each term will be used for specific reference in the following calcu- 
lations. 
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6.1 Evaluation of the first moments of H 

This section shows a detailed calculation of the expected values of the various terms of H in 
Initially one notices that from ( |3.51 ) and ( |4.15 ), with Vt=FpFh, comes 



{FpFh) = 5^5]e[^^^l = W_{l + xp^)\{{l+xhu) = fp{x,y)Mx,y) = f{x,y) , (6.11) 

(1) (2) M 



then, for the first term of Hhf in (6.6) one has 



(12)({/M) (a) {12)(UAI) (a) (1) (a) (p=a,h) 

= EE(^) E e[-i.. = E^'^-^E(^) E (6-13) 

(1) (a) (p,/i) (a) (p-Ma,/i) 

where ip=a) represents all configurations with p particles in which the sp-state (q) is present, 

S = Ua — and t = Ma — niQ, (6.14) 

and the domains (grids) of Ua and Ma are obtained from the domains of U and M by selecting only 
the energies and angular momenta of the {p=a) configurations. When writing the first expression 



of (6.13) the following important relation was used 



E (•••)= E (•••)' (6.15) 

(p-l^a) (P=a) 

which is correct in the sense that the sets of configurations in both sums are in bijective correspon- 
dence with each other (i. e., for each configuration in the sum of the LHS corresponds one and 
only one configuration on the RHS) with the same degeneracies. 

Therefore, the following relation between the degeneracies of (p-l^a) and {p=a) is also valid, 
d{p — 1/a, h, S, t) = d{p — 1/a, h, Ua — Ma — ma) = d{p = a, h, Ua, Ma) , (6.16) 

which express the invariability of a configuration in which a given sp-state (a) is destroyed and 
immediately re-created aftwerwards. The number of configurations and the corresponding degen- 
eracies decrease by the destruction of (a) but they are kept invariable by its subsequent re-creation. 
It is usual to write simply (U-ea) for the energy argument of the degeneracy on the LHS instead of 
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(Ua-ea), but this may lead to inconsistencies when comparing the detailed microscopic grid of the 
various energies. 



Then, using (6.11) results 



■pa 



{a) {p-l,h) 



From (6.15) one can also write 



{p-lf^a,h) 



= EEE^^""""'"'^" E (1) ''-''EE^"*'E^"'"""'"^"^(^^- 1^"'^'^'*) ' (^-is) 

(1) (a) (St) (j=l) (1) (St) (a) 



(a) 



1 + 



(6.17) 



pa 



E E e[^'^^^'^]^e,d(i5 = a,/i,C/„M„) 

(1) (C/„A/„) (a) 



(6.19) 



where the approximation used in the last expression is that, under CAP, the difference in the 
domains of the various 5 and t as a function of (a) can be neglected, then 



and the identity 



EE -EE 

(a) (St) (St) (a) 



E(-)= E (■ 

(St) (C/aM,) 



(6.20) 



(6.21) 



has been used, wihch is a direct consequence of (6.16), again due to the bijective correspondence 



between the domains of {S,t) and {Ua,Ma). This correspondence is not changed by CAP. Expression 
(6.20) is the CAP correspondent of ( |6.13 ) and (6.19) is the correspondent of (6.14). 



In Appendix IV we analyze some additional aspects of these results and some problems that 
may appear if the Laplace transform approach is used in this case. 



For the second term of (6.6) the calculation is analogous and gives 



(/3) (pA) ip,hm 



1 + Xhi3 



(6.22) 



Therefore, 



= EEE^'^""'^^/^ E (1) - Y.Y.-^''^Y.--'''-''''^d{p,h-i^p,s,t) , (6.23) 

(1) {St) (j=l) (1) {St) {13) 



where 



S = U -ep and t = M - 



(6.24) 
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Notice that, under CAP and according to [Appendix IV the sp-state exponential terms in 



Ca and e/s in (6.18) and (6.23), respectively, should be "absorbed" into the integral of the Laplace 
transform if the remaining terms are to be interpreted as the inverse Laplace transforms of the 
corresponding expected values. 



In the residual interaction of Eq.(6.10) the only terms that have non null expected values are 



those indicated by (d), (e) and (/). The calculation is again straightforward, for example for (d) 
one obtains 



(d) = ^ ^ Vai3^s{ph\a+a'^asaj\ph) 



(6.25) 



E 



(6.26) 



(1) (q/357) 



{p-2f^al35^,h) 



(6.27) 



T.y'-' E 

{p,h) {a/3,57) 



{p-2=/=a(35^,h) 



(6.28) 



(p,h) (a/3) 



E 



J5i] 



(p-2=^al3,h) 



(6.29) 



^E^a/3a/g(^pa^p/3) E ^^"^^ ^ E ^ 
(a/3) (p-2,/i) {p-2^aP,h) 



(t^) ^''"^ (l + x,„)(l + x,;,) • 



(6.30) 



(6.31) 



The term corresponding to (e) is 



E E ^"/375(P^I«a^l^/9«7bM 



(-,)(«/3<57) 



(6.32) 
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V'Q,/3-y5(5(p = /i = /3/(5)(5q,,^J/3,5) (6.33) 

(,-)(«/357) 



y; e[^*l (6.34) 



Y.Vapap{xpaXhp)Y.v''"^xP-^ ^ el^*] (6.35) 

(a/3) {p^c,h^P) 

^ (i+x,„)(i+x.,) ' ^^-^^^ 



and the last non null term in ( 6.10[ ) is 



^^)= E Z yo^Pi6{ph\hth+h^bp\ph) (6.37) 

(,-)(a/357) 

X] X] ^°'3T'5'^*^^'^ = "^^'^'^)*^'^".'5'^/9:7 + '^«.7'^/9,5) (6.38) 
(-,)(«/3^7) 



^ K.;3,5(^^)(<^.,.5/3,, + <5.,,5/3,5) e[^*l (6.39) 

(a/357) V ^ / (p,h^a/357) 



= 25^K.;3„/3(x;,„x/,;3)5^y'^-2xP el^*] (6.40) 

(a/3) (p,h) (p,fe7^a/3) 

= 2i;ya/3a/3(x.a../3) + ,,^) + ,,^) > (^-^D 

(a/3) 

If one invokes CAP over the expressions of the expected values and follows a reasoning similar 
to [See. 4 they become the Laplace transforms of the operators' expectated values multiplied by the 



corresponding degeneracies. Consequently, these linear combinations of nuclear degeneracies can 
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be identified with the inverse Laplace transform of the above expressions. This can be achieved by 
simply rearranging the sums intead of doing the expansion of 

For example, from ( 6.27[ ) one has 



(1) ^ ^ ip-2^af3,h) 



and instead of expanding '^^^ write 

(^) = 2EE^"/^"4^) E e[-*U2j;j:e[^'^^^^^]5;K.,.,X:(l)' (6-43) 



(1) M) ^ (p-2^ap,h) (1) (Si) (a/3) (j) 



where d is the degeneracy associated with the energy S and angular momentum t of the configura- 
tions selected by the destruction of sp-states, given by 

S = Uap - e-a- and t = Map - ma - (6.44) 

and Uap and Map are the energy and angular momentum before the destruction of (a) and (/?) of 
the selected configurations in which (a) and (/3) are simultaneously present. Then, 

(d)«2^ E e[^"^*^'^'']^14;3„;3d(p-2/a/3,/i,5,t) . (6.45) 

(1) {U^pM^p) (a/3) 



= 2E E e[^"^*^<^'']^y,^,;3d(p = a/3,/i,C/«;3,M„/3) . (6.46) 

(1) {U^pM^p) (a/3) 

where the identity 

d{p - 27^a/3, /i, 5, t) = d{p = ap, h, U^p, M^p) (6.47) 

was used. 

Note that ( 6.45[ ) can be rewritten as 



{d) ^ Yl yVeI-^^12^14^,^ Yl ^-^^-'d{p-2^aP,h,S,t) , (6.48) 

{pXU) (a/3) {M^p) 

where the sum over U can be approximated as the Laplace transform of the last two sums in the 
CAP limit. 



Analogously, the calculation of (/) yields 
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(1) {UM) {a/3) 



(1) {UM) M) 



(6.49) 



with S and t given by (6.45) and {U ,M) were used as a simplified notation for {Uai3,Ma 



Due to tlie use of CAP in all approximations, the last sums in (6.46), (6.48) and (6.49) 



can be interpreted as the inverse Laplace transform of the corresponding sums over all sp-states 
of the expected values of the residual interaction weighted by the degeneracies for given energy 
and angular momentum. These sums can be used as the definition of the first moments of the 
corresponding operators. [TT] 



6.2 Evaluation of the second moments of H 

The second momenta of the Hamiltonian are also called transition strengths and are con- 
nected with the transition rates in pre-equilibrium models. [U [7] 

Some of these terms do not have a dependence on the excitation energy involving all processes 
with sp-states and therefore are not "on-shell" processes. These "virtual" processes are already 
taken into account by Brueckner's theory and one should be careful with the possible double 
counting of these contributions when calculating the transition strengths. [U [T3] 

Due to the greater complexity involved in these calculations in comparison with the first 
momenta, the detailed calculation of (a|a^) will be presented next in a subsection of its own. 



6.2.1 Terms that increases the number of p and h by 2 



hi Eq.(6.10) the term that increases the number of particles and holes by 2 is (a). The 



corresponding expansion is similar to what was done in the previous section, 

(a|a^) = {ph\a'^a'j^b^b^\ph'J^ph'\b.ybsai3aa\ph) (6.50) 

(-,)(-/3^7) 

= E E \^c.i5-ys\^^{p = = dj)S{p'^a(3)6{h'^6^)6{p\p' + 2)5ih\h' + 2) . (6.51) 
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Now one factor ates the terms in et^*^l that have been selected by the destruction of sp-states, 

to obtain a grand canonical distribution in terms of the energy and angular momentum of the 
intermediary state, 

(«I«^) = E E \Va,,s\' h"'''f''f'A E el'%{p\p' + 2)6ih\h' + 2), (6.52) 

where 

S = U — ea — — es — e-y and t = M — hIq — ni/j — — , (6.53) 

expands 

(a|a+) = J2 y'^'y"'^"' E E + + 2) (6.54) 



and, similarly to ( |6.18 ), changes the sums over (p=a/3) and {h=S"f) into sums over {p-2j^af3) and 



{h-2^6'y), as suggested by the new distribution in terms of S and t, 

= E \ya^-ys\\xpaXp^XH5Xhy) ^ /-'y'^'x^-'x^'' ^ e[^%(p - 2|p')<5(/i " 2|/i') (6.55) 

^ ' \p',h') \ p'j^cl3,h'^S-f J 



_ ^ 2/ Uf.[^ + Xpt^x'j^,^)Uu{'^+XhuX'f^,^) 

~ / J l^alS-yS] {XpaXpi3XhSXh^)— ■ - r— ■ -. — - ■ - r— ■ - r . (6.57) 

As we have seen before, this expression can be considered, under CAP, as the Laplace transform 
of the expected value of the residual interaction times the corresponding degeneracy, for each 
(a,/3,7,(5). 

To obtain the expression of the corresponding inverse Laplace transform one may procede as 



in (|6.43|), then from (|6.52|) comes 

XpaXpfljXfigX}i-y \ 



{aUpib:^\b,bsapa^) =E E l^./^7^n 2 2 ' E ^^'MpIp +mh\h' + 2) (6.58) 
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dd' 



(1) {,^,1,) (»l3Sy) {j,k) 



(6.59) 



where S'=C/'=C/-(ea+e/3+e5+e^) , t=M'=M- (mo +m/3 +1115 +m^) and 

S' = U' + ea + + es + £7 and t' = M' + iiIq, + + + . 



(6.60) 



Here {U,M) are the energy and momentum of the selected initial configuration in which the 
sp-states (a,/3,7,5) are present and (U' ,M') are the energy and momentum corresponding to the 
selected intermediary configuration in which the same sp-states are supposed to not be present. 
Therefore, the domains of {U,M) and {U',M') have the "same number of points"]^ and correspond 
to a simple displacement of each other. 

The identity between the energies, {U,U') and {S',S), and the angular momenta, {M,M') 
and {t',t), is a consequence of the definition of (a|a^) if one is not considering off-shell transitions. 

Due to the supposed orthonormality of the single particle basis, for each configuration 
\p — 2,h — 2) associated with the energy S and momentum t, one cannot have more than one 
\p'h') satisfying the condition <p — 2,h — 2\p'h'>=l. Therefore, the last sum over (j = l,d') is 
equal to 1 and the entire sum over the intermediary states \p'h') reduces to only one term. 

On the other hand, the sum over (k = l,d) yields the nuclear degeneracy associated with 
{p,h,S,t) and the expression of (o|a^) can be rewritten as 



(a|a^ 



^^el^^^l E \VaP'y5\^dip-2^a(3,h-2^6^,S,t) , 



(1) (St) (a/357) 



(6.61) 



where now 



AU,M] 



J-/3(7-7Af-/3l/'-7M'] 



CAP 



>2[-/3C/-7M] 



(6.62) 



due to energy conservation, i. e., the probability distribution after the transition, in which the final 
nuclear state equal to the initial one, is the square of the grand canonical probability of existence 



of the initial state. The factor of 2 is non important in the CAP limit and (6.61 ) becomes the usual 



Laplace transform expression as we saw in Sec. 4 (equations (4.11) and (4.14)). 



Equation (6.61) is analogous to the first moment expressions in (6.47) and (6.50). There are 



no essential differences in both cases because the present analysis reduces to on-shell processes in 
which the intermediary states do not play an important role on the calculated expected values. 



^This is correct whether the domains are continuous or discrete grids. 
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In Appendix V one shows that (6.61) can be written as the foUowing convolution involving 



densities of the nuclear states associated with the given degeneracy, 

(a|a+) ^ J^^el^*^]^ j deV{2,2,€,:a) u}{p - 2,h - 2,U - e, M 



(1) (St) (m) 

where e=ea+ei3+e-y+es, m=nicj+ni^+m^+in5 and 



y(2,2,e,m)= J dx dy de4\Vai3'ys\'^ oj {1,0, e — x,m — m:^) 



(nil, my ,1115) 



X a;(l,0,x - y,mx - my)a;(0, l,y - €4, my -ms)uj{0, 1,64,1115) 



(6.63) 



(6.64) 



with y=es+e4, x=e2+y, mx=m^+niy and niy=ni^+ma. 



6.2.2 Terms that increases the number of p and h by one. 

The strengths of the transitions that increase the number of both p and h by one correspond 



to the terms (6) and (c) in Eq.(6.10) and the corresponding terms that decrease p and h by one are 
(g) and (/i), respectively. 

The transition strength corresponding to term (6) can be calculated in a totally analogous 
fashion as for (a|a^) in the previous section, 

(6|6t) = ^ ^ \Vap^s?{ph\a^a^hja^\p'h'\p'h'\a^hsaisaa\ph) (6.65) 
= E E l^"/375l'^b = oiP^l)5{h = 5)5{p' = ^^aP)6{h'^5)6{p\p' + l)5{h\h' + 1) . (6.66) 

(-,)M67) 

Now one factorates the terms in e^^*^] that have been selected by the destruction of sp-states both 
in {ph) and {p'h'), redefining the grand canonical distribution in terms of the energy and angular 
momentum of the intermediary state, 

U' = S = {U — ta — — + 67) and M' = t = {M — hIq — — + m^) , (6.67) 



then, 



(fe|^^) = EEl^-/^.^l V'"X"'^i E eiSMp-2\p'-mh-l\hr (6.68) 

(1) (a/357) V ^ Mr"5^;i'^r') 
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Therefore, {U,M) correspond to the selected initial configuration in which the sp-states {a, (3,6) are 
present and (7) is not, and {U' ,M') to the selected intermediary state in which the sp-state (7) is 
present and {a, (3,5) are not. 

Notice that the term (Xp,^/ x')=e^^^''^^^'' is not present in e'^*^] because (7) is not supposed 
to belong to the initial configuration (ph) , but the conservation laws imply that and must be 
added to {U ,M) to produce {U',M'). Then, to write the grand canonical distribution in terms of 
{U',M') one must multiply the entire expression by [x'^i^/x'). 

Now one expands obtaining 

m = Y^ yWx'^X \V^,,s? h''"'f'fA Y: el'Mp-2\p'-mh-l\h'){6m) 

\p',h'J ^ ' \p'=-,T^al3,h'j^s) 



\p',h') \ p'-l^a/St.h'ytS ) 

and rewrite some terms a using the deltas 

= E \Vap,s\'{x,^x,^x,sx',,,) E {yy')^''-'\xx')(^-'^ Y 



, ,2, , , UJ'^ + Xp^,x'p,JUA'^+^hux',,,^} 

y , IValS-ySl [XpaXppXhSXpi^) j r (6.71) 

(a/357) + XpaX'p^) { 1 + X^px'^^ J (l + Xh&x'^^) (l + Xp^x'^^) 



To obtain the expression corresponding to the inverse Laplace-transform one procedes as in 



Eq.(6.58), 



{ata^pb^a,\apsapa^) = Y E l^./^^^P r'^^'r''"^"'" E -^''HpW + mh\h' + ^.72) 



dd' 

Ew+i)^(^/^'+i) E E \yo.Ms?^(u - s')6{M - t')E(i) • (6-73) 

(1) (s5*,) (°/3'57) (i.fc) 
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where 

S' = [U' + + + — e^) and t' = (M' + iiIq, + + 1115 — m^) . (6.74) 

As happened for (a|a^), the identificatfon of the energies, {U ,U')={S' ,S), and the angular 
momenta, {M,M')={t',t), is a consequence of the definition of {b\b^) where only on-shell processes 
occur, therefore, the sum over intermediary states (p'h') reduces to a single factor of one due to 
the various constraints relating them to the initial and final states. Then, 

= ^^e[^*^l \V^p^s\^d{p-2 + l^a/3j,h-1^6,S,t) . (6.75) 

(1) (5,i) (a/357) 

Here d{p — 2 + l^af3^, h — lj^6, S, t) is the degeneracy associated with the configurations 
obtained by the destruction of 2 "p" sp-states followed by the creation of 1 "p" sp-state, plus the 
destruction of 1 "h" sp-state and ( 6.75[ ) could also be written, more explicitly, as 



(6|6+) = ^^el^*^] ^ \Vo,^^s\^dip-a-f3 + ^,h-5,S,t) . (6.76) 

(1) {S,t) (a/357) 

One notices that "p"s and "/i"s are created and destroyed independently and, although 
d{p — 2 + l,h — 1, S,t) is in general different of d{p + 1 — 2,h — 1, S,t) as functions of S and t, the 
sums of the degneracies over all configurations are the same, 

NOC{p - 2 + 1, - 1) = Yd{p -2 + l,h-l,S,t) (6.77) 

(5,t) 

= Yd{p +l-2,h-l,S,t) = NOC{p + 1 - 2, /i - 1) , (6.78) 

(5,t) 

where NOC stands for "number of configurations". 

Therefore, under CAP it becomes reasonable to drop the explicit dependence on the sp-states 



in (6.75) and (6.76) and rewrite them as 

(6|6+)^^^j;^e[^^^] \V.p,8?d{p-2 + l,h-l,S,t), (6.79) 

(1) (S,t) (a/357) 

where the contributions of terms with a given {S,t) are the same even if they correspond to different 
sets of sp-states. 



In Appendix V one shows that (6.79) can be rewritten as a convolution 



|5+) ^et^*^] Y \VaP'y5?{d{p-a- l3,h-6,S,t)-d{p-a- ^--i,h-5,S,t)) (6.80) 
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^e[^^^^l[^ / d ey^(2,l,a;,mx)w(p-2,/i- l,C/-x,M-mJ 
/ deV{2, 2, e, m) u{p - 2,h - 2,U - e, M - m)] , 



(6.81) 



where x=(ei+e2+e4), y=(e2+e4), nix =(mo +111^+1115) and niy =(111^+1115), 

T4(2, l,x,mx) = / ^"(^^5) w(l,0,x - y,nix -niy) w(0, l,y - e4,niy - ni5)a;(0, 1, €4,1115), (6.82) 



and 



= j desy^JVaiBysl'^ UJ{1, 0, £3,111^) , 



y(2,2,e,m)= J dx dy de4,\Va/3'ys\'^ ^{^,0,e — x,ni — 

(nix jUiy ,m j ) 

X a;(l,0,x - y,m^ - my)a;(0, l,y - e4,niy -ms)u:{0, 1,64,1115) 



(6.83) 



(6.84) 



with e=(ei+e2+e3+e4), in=(ma +111/3+111^+1115), x=(e2+e3+e4), mx=(m/3+ni^+ma), 7/=(e3+e4) and 
my=(ni^+mQ). 



6.3 Terms with off-shell contributions 



The various terms in (6.30) are defined as sums over sp-states and in some cases both 



Hermitean conjugates of a given sp-operator, Upa and Opa, act on the same configuration and the 
corresponding sp-energies become absent from the final accounting of the energy conservation for 
the total transition. In these cases the process of creation and destruction of these excitons is 
off-shell with respect to the entire transition. 



For example, in the case of the expression of {b\b^) in (6.69) with 7 replaced by a in (6) and 
(6+) one has 

(&^^a|6+_^„) = ^ 'Yl I Vai3s \'^{ph\a'^a'^bjaa\ p h' \p' h' \ a+ 65 a pap, \ ph) (6.85) 
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which gives the following constraints over configurations 

= \^o.0s\'^^ip = aP)S{h = 6)6{p' = a^p)5{h'^6)6{p\p' + l)6{h\h' + 1) . (6.86) 

Now one factorates the terms in e^^^^ selected by the destruction of sp-states, 

ib.^a\bUa) = E E hT'f'^ ) E -''Mp - 2\P' - mh - l\h') , (6.87) 

where 

S = {U -ei3- es) and t = (M - - 1115) , (6.88) 

and, therefore, the sp-cncrgy and momentum {ea,^a) are not present in the formal conditions of 
energy and momentum conservations. 

Now one expands rewrite the terms using the deltas, 

= E y'-'y"'-" E i^^/^^i' hT'f'" ] E ^'''^(p - ^\p' - mh - m (6.89) 

/ p,h \ (aBS) \ y / ( P=oil3,h=5 \ 

\p',h') ^ ' \p'=a^l3,h'^6) 



E \yapsWxp^XppXusx'j,,^)Y, y''-'xP-^y"''x'P'-' E ^^''^^(P - - mh - l\h'),i6.90) 

(a8S) ( P'f^ \ (p-2^afi,h-l^S\ 

\p' ,h' ) \ pi -l^aP,hl^S ) 



to obtain the Laplace transform of the expected values of these operators, 

= El^«/35l'(wp«4'a)E(^^')^""'Hyy')^'"'^ E ' 



[St] 



= 2^\Val3&\ [XpaXppXhSXp,^)- -y — (6.91) 

(a/35) (1 + XpaX'pa) [1 + Xp^X^^ J (1 + XhSXf^g) 

As before, the expression for the inverse Laplace-transform can be obtained by rearranging 
the terms and using CAP, 
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dd' 

CAP 



where 

5' = {U' + efs + es) and t' = {W + m^j + m^) , (6.94) 

which reduces to 

« EE^^'"'"' E - 2 + l/<5,5,t) , (6.95) 

(1) (5,t) (a/3<5) 

and can be written more explicitly as 

(6^^,16+^,) « EE^""*'' E \^o.ps\^dip -a-/3 + a,h-6,S,t) , (6.96) 

(1) {S,t) (a/35) 



In Appendix V| one shows that ( |6.96[ ) can be rewritten as a convolution 



(6^^,|6+^„) « ^^e[^^^] y"dx^y,(2,l,x,niJw(p-2,/i-l,;7-x,M-mJ (6.97) 



(1) K) 



where 



K(2, l,x,mx) = J dydes ^ |V;^/3^5ptj(l, 0, x - i/,mx - niy) 



(%"1<5) 

xw(l,0,y - e3,niy - ms)uj{0, 1, €3,1115) , (6.98) 



and 



X = ea + ^13 + es and nix = hIq, + + 1115 . (6.99) 



When acting over their corresponding Hermitean conjugates the sequences of sp operators 
do not need to be the exact conjugates of each other as it has happened so far, for example, the 
term obtained from the previous one with a replaced by 77^0 in (6+) is defined as 

{by-^a\ba^^)^^^= ^ ^ \Vai3^sf{ph\aialblaa\p'h''^ph'\a\b5a[say\ph), (6.100) 

(-,)("/357) 
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and it can be treated similarly. 

Notice that the selection of sp-states must correspond to the action of the Fock operators 
acting both on the nuclear "ket" and "bra" states, for example, in the present case the operation 
of (aao'pbgaa) on {ph\ produces the selection of the sp-state a which is not present due to the 
operation of (a^^^^a^a^) on \ph) , etc.. This is a simple rule that results from the conservation laws, 
because the intermediary state is unique. 

Then, the selections are represented by the following constraints over configurations, 
= E E \Vaf)js\''Sip = jPa)d{h = 5)5{p' = a^^/3)S{h'^5)dip-2\p' - l)d{h-l\h') .(6.101) 

Now one factorates the terms in e^'^^ selected by the destruction of sp-states 

(1) (a/357) V ^ /( r^^"^pD 

expand and rewrite some terms a using the deltas, 

= Kp-ysl^ixpjXp^XpaXhsXp^^x'p,^) Yl {xx')^-^\yy'i^-^'^ Y 

{Xp^XppXpaXhSx' i^x' , )W ( 1 + Xpi^x' I \ (l+XhuX'f^,^) 

= J2 l^«/37<5l 7 — — ■ (6-103) 

(a/357) (1 + ^PcX'pa) (l + Xp^x'^^j (l + XhSx'f^) (l + Xp^x'p^) 



The expression corresponding to the inverse Laplace-transform is obtained as before, 

("I"*),.,. = E E ( ^-'''"■yjf''--''-' ) - 3b' - mh - i|/.'i6.io4) 

(1) (a/357) V / (P-3^'(l3o,,h-l:^S\ 



dd' 



(1) («/357) m 



which reduces to 



(6|6+)^^, « ^J^el^^l Y d{p-S + 2^p,h-l^S,S,t), (6.106) 

(1) (5,t) (a/357) 
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and can also be written as 

« 5] J^el^^'^l dip-a-j-P + a + ^,h-5,S,t) , (6.107) 

(1) iS,t) (a/357) 

where now both a and 7 are canceled out in the expressions of energy and angular momentum 
conservation, 



S={U-ep- es) 



and 



t = (M - - m^) 



(6.108) 



and 



S' = [V + ep + e&) and t' = [M' +111/3+ m^) 



5.109) 



Therefore, the transitions to and from sp-states a and 7 are described as off-shell. 



The degneracies in 6.75 and 6.95 should have similar magnitudes as functions of U and M, 



but they sually would be considerably larger than (6.106), despite the fact that all them correspond 



to configurations with (p-l,/i-l) excitons, which should have a direct reflex on the magnitudes of 
the corresponding transition strengths. 



In Appendix V one shows that (^6.106) can be rewritten as the following convolution, 



where 



(6^^„|6+^^)^_,^ « /'a!eJ^K(3,l,e,m)^^(p-3,/i-l,C/-e,M-m) (6.110) 

(m) 

V;(3, l,e,m) = / dxdydei ^ |V'Q,/3^5pa;(l, 0, e - x, m - m^) 



xa;(l,0,x - y,mx - my)a;(l,0,y - £4, my - 1115) w(0, 1, £4,1115) , 



.111) 



e = + e/3 + + £5 and m = iiIq + + + 1115 



X = + + and nix = + + 



5.112) 
5.113) 



and 



y = + and niy = + 1115 , 



(6.114) 
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6.4 Other terms with on-shell transitions 



At this point it is clear that all transitions can be treated in essentially the same way to 
obtain the final expressions in term of the degneracies, with or without CAP. 

The basic steps are the fatorization of the exponential terms associated with the sp-states 
destroyed in the transition and the rearrangement of the sums using the deltas to express them 
in terms of the quantities corresponding to the intermediary states. We will use this rule in the 
following analysis. 



The term corresponding to (c|c+) in (6.10) is 



{c\c+)=Y, (P^\''^c.blb\bp\p'h'lp'h'\blb^bsa^\ph) (6.115) 

which imply the following constraints for configurations, 

(p = a), (/i = <57//3), {h' = /3/57), (p'^a) (6.116) 

and the following terms factorized from e^^^^l, 

{xpaXhSXh-yXj^i^) . (6.117) 

Then, using CAP the following final expression is obtained 

(c|c+) « J]^e[^*^l ^ \V^^^s\^d{p-l^a,h-2 + l^p6j,S,t) , (6.118) 

(1) {S,t) (a/357) 

or more explicitly, 

(c|c+) « J^el^A^l \V^p^s\^d{p-a,h-6-j + /3,S,t), (6.119) 

(1) {S,t) (afiS^) 

where the interemediary energy and momentum are 

U' = S = {U - ea- es - + 6,3) and M' = t = (M - - 1115 - + m^) , (6.120) 



and one can show that (|6.120) can be rewritten as the following convolution (see Appendix V) 



(c|c+) « Y^iUM] \v^^^^\^[d{p-a,h-6--f,Si,t)-dip-a,h-(3-j-6,S2,t)] (6.121) 
is,) 
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^e[^^l[J] fd eV^{l,2,x,m^)uj{p-l,h-2,U-x,M-m^) , 

j (iey(l,3,e,m)a;(p- l,/i-3,f/-e,M-m)] , (6.122) 



where 



(m) 



Si = {U — Ea — ^■y — €5) and ti = (M — niQ, — — m^) , (6.123) 



S2 = {U — ea — — — es) and t2 = (M — hIq, — — nLy — m^) , (6.124) 



a;= (61+63+64), y=(e3+64), nix =(mQ, +111^+1115) and my =(111^+1115), 
Vx{l,2,x,m^)=Y / dyd64 V'(Q,^5)a;(l,0,x -?/,mx - my) a;(0,l,y- 64, my -m5)(x;(0, 1,64,1115), (6. 125) 

(my, ma) 

V(ajS) = / de2Y\Va(3jsf UJ {0,1, 62, mp) = Y\Vaf}^sf , (6.126) 

and 

y(l,3,6,m)= ^ Jdxdyde4^\Vaj3^s\'^u;{l,0,e-x,m-mj,) 

(mj,my,mj) 

X a;(0, l,a; - y,mx -my)a;(0, l,y - 64, my - m^) w(0, l,64,m5) , (6.127) 



with e=(ei+62+63+64), m=(mQ,+m^+iiLy+m5), ^=(62+63+64), mx=(m^+m^+mQ,), j/=(63+64) and 
my=(m^+mQ,). 

The terms corresponding to {c\c^) 0^5 with off-shell transitions for (/3) and (S) in c and c''' 

are 

{c\c^)p^S=Yl J2 ^P''\''ibiKMph'lp'h'\b+b^bsac,\ph) (6.128) 
(,-)M5t) 

which imply the following constraints, 

(p = a), {h = 5jP), {h' = PS^J), (p'^a) (6.129) 



and the terms factorized from e^^'^^ are 

(^XpaXjiSXfi^Xlif^Xf^^Xf^pj , 

giving the foUowing expression in terms of the degeneracies of the intermediary states, 

(1) (5,t) (a/3<57) 



or 



Y^Y^e^UM] J2 \v^^^,fd{p-a,h-f3-6-^ + (3 + 6,S,t) , 

(1) {S,t) (a/357) 

where 

U' = S={U-ea- e-y) and M' = t = {M - - e^) . 
Finahy, equation ( |6.133[ ) can be rewritten as. 



{(l),S,t) (m) 



where 



V;(l,3,e,m) = J dxdydei ^ |V'a/3^5|^a;(l, 0, e - x, m - nix) 



/rnxmy \ 



xa;(0, l,x - y,mx - my)a;(0, l,y - e4,my - 1115)^(0, 1,64,1115) , 
e = + 6/3 + + 65 and m = iiIq, + + + 1115 , 

X = ep + e-y + and nix = m^j + + 



and 



y = 6^ + 65 and niy = + ni5 . 
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7. Results and conclusions 



We have shown in Appendix IV (see also related notes in App.I.l) that if one relies on the 
Laplace transform as the main formalism for the transition strengths (TS), some microscopic tran- 
sitions are not correctly described, or not well defined. Therefore, a direct microscopic formalism 
is preferred. 

The formalism presented in this paper is a consequence of the microscopic definitions of the 
first and second momenta of the nuclear Hamiltonian with a residual interaction term. It provides 
a direct and intuitive understanding of the TS (for given number of excitons and energy of nuclear 
state) and its dependence on the strengths of less complex states, which can be expressed in terms 
of convolutions of these functions, in a similar way to the traditional relations obtained for nuclear 
densities. These expressions have been used for the evaluation of the TS and the densities and 
compared with the corresponding parameters of the exciton model (EXM). 

The Model Space was defined by arbitrarily fixing the Fermi level, Ep, as greater than or 
equal to the minimum energy necessary to accomodate all protons and neutrons of the system in 
the fundamental state, and by considering the levels of the basis above and below Ep, to describe 
the fields of holes and particles respectively, up to a maximum energy per exciton, -Emax- The latter 
was then fixed arbitrarily as a model parameter. 

Various "sizes" of the Model Space were considered by varying E'^ax and Ep. For sake of 
simplicity a not too large nuclide was used, ^"Ca, to permit to test the model with relatively small 
runtime computation, but even so the number of configurations for increasing n increased very 
fast making it unpractical to perform direct calculations for n greater than 4. Calculations with 
larger number of excitons, would probably need an approximate statistical approach and should be 
attempted only for larger nuclei, as the consideration of large n in small nuclei is usually of little 
physical interest, or otherwise if a specific practical situation demands it. 

For example, for ^"Ca, "n=3" and using a single "personal computer" (PC) the number of 
configurations is of the order of 10^-10^ with computation time of a few minutes up to a couple 
of hours, depending on the considered transition and the size of the Model Space. For "n=4" this 
number approximately doubles and the maximum runtime reaches many hours. For "n=5" the 
number of configurations can be 10 times greater and the estimated runtime would reach many 
days. In addition, the calculated strengths increase by two or more orders of magnitude for an 
increase of An=+1 making it more difficult the direct computation of the resulting large numbers 
while keeping the same numerical precision. Nonetheless, the results seem clear and meaningful 
enough even with 4 excitons or less and no more than that will be considered in the present analysis. 
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The idea of convolution numerically results from the fact that all strengths, for all numbers 
of excitons, energies and momenta can be obtained from a core of elements of the transition matrix 
with no constraint over the microscopic initial state. This permits to obtain a basic set of TS 
for the simplest configurations for which the transition is possible and the TS for states of higher 
excitations by doing the convolution of the basic TS with the adequate nuclear densities. 

The core of matrix elements can be calculated independently for given maximum excitation 
t^iax, given Ep and the assumptions about the single particle level density, then its components 
are used to obtain the transition strengths systematically for different number of excitons, n, as 
a function of U . In addition, the initial configuration is constrained by the necessary presence of 
the sp-states that will be "destroyed" during the collision process and the final by the exclusion 
principle for the newly created sp-states. 



For example, in Eq.(6.122) it was shown that the sum of TS designated by (c|c"'") can be 
written as 

(c|c+)= ^ {ph\aib\h\hp\p'h'lp'h'\h\h^ha^\ph) , 

~ ^e[^^^][^ / dxT4(l,2,x,mx)a;(p- l,/i-2,C/-x,M-mx) , 

/ deV{l, 3,e,m)uj{p-l,h-3,U -e,M - m)] , (7.1) 

W 

where V^(l,2,x,mx) and 1^(1, 3, e,m) are the TS corresponding to initial configurations with 3 excitons, 
{p=l,h=2), and 4 excitons, (p=l,/i=3), respectively, which are the simplest possible configurations 
in this case. The strengths V^(l,2,2;,mx) and V{l,3,e,m) are defined as sums over the elements of 
the core of transition matrix elements that are independent of p, h, U and M, therefore, they are 
functions of the Model Space but not of the level of complexity of the excited nuclear state. 

The densities that enter in the definition of the transition strengths are calculated indepen- 
dently either by combinatorial analysis, for the initial nuclear configurations, or by direct counting 
of the available states in connection with each valid transition. 

Using this prescription, a computer code called TRANSNU[8j was created to calculate the 
TS. We also assumed a "billiard balls" model for the microscopic description of the system of 
excitons, i. e., a quasi- free movement of the excitons in the nuclear matter and interaction by 
a macroscopic-type "collision", i. e., a distance independent satationary redefinition of the local 
phase space of the interacting sp-states with conservation of energy and angular momentum. 
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Therefore, the cohision was described by the simplest possible meachanism, which is the 
free propagation of the of the initial and final exciton states in the Fock space, followed by the 
destruction of part of them and the creation of the final ones constrained by the conservation laws 
and the exclusion principle. No assumption regarding the dependence of the interaction on the 
relative distance of the excitons during the collision was made, therefore, one may also regard this 
as a type of "black box" description for the microscopic two-body interaction. |15j 

The Greeen function for the pair of "colliding" excitons was the propagator of the non 
interacting two-body Hamiltonian|16j 



e 



ik\fl~fi\ 



kY,{2l + l)2i{kr<)h\^\kr>)Pi{cosa) (7.2) 



ri — ro 

{1=0) 



and the sp-state basis was the set of Harmonic Oscillator (HO) wave-functions, following the pre- 



scription of Brussaard and Glaudemans[5] and as wave number for the stationary excitons in (7.2) 
we used the inverse of the characteristic length of the H0|5] 

k = . (7.3) 

The transition matrix elements were then obtained by direct integration over the entire space of 
relative positions of the colliding pair of excitons. 



7.1 Numerical results 

Inasmuch as the HO basis can be considered as a non phenomenological assumption, the 
above definition for the propagation of the colliding excitons can be considered as non phenomeno- 
logical. With this prescription, the integrals of the transition matrix often produced very large 
numbers due to the non normalization of the two-body propagator. 

Then, taking into account the intrinsically hypothetical nature of the HO description, an 
ad hoc model parameter was introduced to obtain always finite matrix elements that could be 
compared among themselves as functions of the nuclear excitation, as well as with the analogous 
parameters of other models. In this sense the results we present now have only qualitative validity, 
but they can be made more physically meaningful by fitting of the adquate experimental data. 



The various transitions of Eq.(6.10) were classified according to their An and those with 



negative An were neglected, in agreement with the usual procedure of semi-classical analyses. The 



transitions in which the number of excitons increase correspond to the terms (5) and (c) of (6.10) 
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which have the same TS as their Hermitean conjugates, terms (g) and (h) respectively. In the 
numerical analysis these terms have been divided into the following subprocesses 



case 


representation 


description 








1 


(2100++1000) 


proton-h-propagates 


and 


proton-ph-pair 


is 


created 


2 


(1011++1000) 


proton-h-propagates 


and 


neutron-ph-pair 


is 


created 


3 


(1110++0010) 


neutron-h-propagates 


and 


proton-ph-pair 


is 


created 


4 


(0021++0010) 


neutron-h-propagates 


and 


neutron-ph-pair 


is 


created 


5 


(1200++0100) 


proton-p-propagates 


and 


proton-ph-pair 


is 


created 


6 


(0111++0100) 


proton-p-propagates 


and 


neutron-ph-pair 


is 


created 


7 


(1101++0001) 


neutron-p-propagates 


and 


proton-ph-pair 


is 


created 


8 


(0012++0001) 


neutron-p-propagates 


and 


neutron-ph-pair 


is 


created 



and the terms which keep the number of excitons constant in Eq.(6.10), i. e. (d), (e) and (/), were 
divided into the following subprocesses: 



case 


representation 






description 


9 


(1100++0011) 


neutron-ph- 


-pair is destroyed and proton-ph-pair is created 


10 


(2000++2000) 


scattering 


of 


two proton-h excitons 


11 


(1010++1010) 


scattering 


of 


a neutron-h and a proton-h 


12 


(0020++0020) 


scattering 


of 


two neutron-h 


13 


(0200++0200) 


scattering 


of 


two proton-p 


14 


(0101++0101) 


scattering 


of 


a neutron-p and a proton-p 


15 


(0002++0002) 


scattering 


of 


two neutron-p 


16 


(1100++1100) 


scattering 


of 


a proton-ph pair 


17 


(1001++1001) 


scattering 


of 


a proton-h-neutron-p pair 


18 


(0110++0110) 


scattering 


of 


a neutron-h-proton-p pair 


19 


(0011++0011) 


scattering 


of 


a neutron-ph pair 



In the following, the results of cases 1 to 8 are designated by VPLl results. The 9th case 
also corresponds to a pair creation process and is designated isolately by VPL2. At last, the cases 
from 10 to 19 are collectively designated by VOO. In the graphs shown below the contributions of 
VPL2 and VOO were added and compared with VPLl. 

One may question the physical meaningfulness of the comparison of the microscopic calcu- 
lations with the EXM, having in sight the important criticisms of the model made by Pompeia 
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and Carlson. [17] In fact, many of the features of the EXM result from statistical approximations 
that can be more or less physically precise, specially regarding the hypothesis of the attainment 
of equal occupation probabilities for all configurations at each stage of the increasing complexity 
chain, i. e. the attainment of perfect configuration mixing in the ensemble of nuclear microstates 
or "equilibrium hypothesis" (EH). 

One of the arguments against EH is that the transition rates (TR) of the exciton model ("Aq" 
for transitions that keep the number of excitons invariable and "A+" for transitions that increase 
it by 2 units) indicate that "Aq", as predicted by the EXM, is not large enough, in comparison 
with "A+", to warrant the achievement of equilibrium between any two steps of the chain. Using 
a Monte Carlo calculation it was determined that Aq should be at least three orders of magnitude 
greater than A+ to validate EH.[T7] 

On the other hand, the EXM and the hybrid model are two largely used tools in nuclear 
data evaluation due to their simplicity in terms of easiness of implementation and allowance for 
improvement in comparison with the observed data. Then, one may analyse the important param- 
eters of the microscopic formalism by direct comparison with well established EXM results, like the 
TNG model code[l2j, to have an idea of the quality of the estimates of the present formulation. 

The microscopic calculations using TRANSNU, for different exciton numbers and assump- 
tions about the energy level structure of the sp-states (SPL), clearly show that for all exciton 
numbers and SPL densities the TS corresponding to the "A+" of the EXM, "t+", can in fact be 
more than 3 orders of magnitude greater than "to" (the analogous of "Aq" of the EXM) in the 
region of not very small and not very large excitations. Then, one of the basic arguments against 
the EH would be invalidated by the present microscopic calculations. 
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Transition Rates (X) and Nuclear Densities (J) 
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Figure 1. Transition strengths and nuclear densities as functions of the excitation energy for 
transitions that keep the number of excitons constant, for an initial configuration with n=3. 
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Figure 2. Transition strengths and nuclear densities as functions of the excitation energy for 
transitions that increase the number of excitons by 2, for an initial configuration with 1 exciton. 

A typical example of this relation is shown in figures [T] and [2| for to and t+ respectively, for 
configurations with n=3. 

In this case to and t+ have variable ratios as a function of the nuclear excitation, E^x, but 
in the region of maxima to is more than 4 orders of magnitude greater than t+, which according to 
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Ref.[T7] would be enough to validate EH. 



Note that the maximum of the densities in Fig. 2 is only 2 orders of magnitude smaller than 



the maximum of Fig. 1 , indicating that the details of the collision process are the main causes of the 
differences between to and t+. For n=4, to ^-nd are shown in figures |3] and |4| respectively, and 
one notices that the ratio between and t+ is smaller but still greater than 3 orders of magnitude 
in the region of maxima. 
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Figure 3. Transition strengths and nuclear densities as functions of the excitation energy for 
transitions that keep the number of excitons constant, for an initial configuration with 4 excitons. 

Notice that the agreement with EXM is restricted to a region of not too high or too low 
excitations and becomes worse for higher n, suggesting that the EXM should be revised in these 
cases. Our calculations indicate that the points at which the TS fastly decrease to zero depend on 
the level structure of the Model Space and n. 

The details of the factors involving the angular momentum sums have been given in a previous 
work Ref.[H] and the present calculations show that, although the ratio between to and t+ is smaller 
in the region of maxima, when these factors are fully considered, it is still greater than 3 orders 
of magnitude for all exciton numbers and sizes of the Model Space considered in this study. This 
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Transition Rates (X) and Nuclear Densities (J) 
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Figure 4. Transition strengths and nuclear densities as functions of the excitation energy for 
transitions that increase the number of excitons by two, for an initial configuration with 2 excitons. 

indicates that these factors do not play an important role in the definition of the TR and could be 
replaced, with good approximation, by fixed phenomenological parameters. 

In the microscopic results both the transition strengths and the densities show strong oscilla- 
tions from one energy to the next, reflecting the existence of regions with very different degeneracies 
even for relatively close excitations. These differences have essentially a combinatorial nature and 
the physically meaningful degeneracies are those associated with the local maxima of the degeneracy 
curve, corresponding to the most probable configurations. 

For example, [Fig. 5 shows these oscillations for the same transition elements shown in Fig. 3 



The results of |Fig. 3 were obtained by fitting a smooth curve over the rapidly oscillating microscopic 
results, by neglecting the local minima and constraining the area below both curves to be the same. 
Therefore, the smoothed curve in Fig. 5[ although approximate, follows closely the local maxima 



of the rapidly oscillating one and is also a physically meaningful description. 

In summary, the main results and conclusions of this work are, first, to show a precise direct 
microscopic formalism (DMF) connecting the definition of the momenta of the PE Hamiltonian 
and the usual TR parameters of the semi- classical models (EXM and hybrid models), revealing the 
connection with the traditional approach based on the Laplace transform and to notify about some 
specific cases in which the latter gives physically incorrect results. 

Second, the TS calculated with DMF show good comparison with the corresponding semi- 
classical parameters in a certain interval of not very low and not very large excitations, becoming 
smaller than the semi-classical estimate outside this interval. For increasing n the relative decrease 
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of the TS in comparison with the semi-classical TR is faster. We interpret these results as indicating 
that the EXM tends to describe better PE states with low excitation because of the competition, 
for higher excitations, with other processes like PE emission and the possibility of transition to the 
compound stage. 
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Figure 5. The same as Fig. 3, showing also the rapidly oscillating results of the microscopic 
calculations and the smoother curve that represents them. 



A third conclusion is that in the region of good similarity with EXM the ratio between to 
and t+ is usually greater than the minimum that permits to validade the "equilibrium hypotesis" , 
as analyzed in Ref.[T7], as one of the central aspects of the EXM. This indicates that the estimate 
for the interaction parameter |Mp of semi-classical rates of transitions should be revised. 

Fourth, the details of the angular momentum coupling play no major role in the above 
results and yield closely the same strengths and ratios between and to ^-s those obtained with 
an adequately defined constant. 

At last, the results shown in the previous figures correspond to a modified HO basis in which 
an arbitrary constant spacing between sp-levels has been assumed. Without this assumption the 
comparison with EXM becomes not so good as we see in Figj6]for transitions with An=0 and FigjT] 
for transitions with An=+2. 

Nonetheless, despite the remarkable difference with usual EXM estimates graphs [6] and [7] still 
show a ratio between to and t+ that is still compatible with EH in the region of maxima. 
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Transition Rates (X) and Nuclear Densities (J) 
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Figure 6. Transition strengths and nuclear densities as functions of the excitation energy for 
transitions that keep the number of excitons constant. It includes the rapidly oscillating results of 
the microscopic calculations and the smoother "linearized" curve that represents them, using the 
HO sp-states levels for n=4. 
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Figure 7. The same as Fig. 6 for transitions that increase the number of excitons by two. 
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Appendix I - A few comments on the Darwin-Fowler method 



Using the Darwin-Fowler formalism the level density can be defined as the pole of the grand 
canonical generating fmiction f{x,y) divided by x^~^^y^^^ [31 H] and, consequently, it can also 
can be interpreted as the inverse Laplace transform of a linear combination of nuclear densities for 
given number of excitons (A) and excitation energy (?7). Here we present a detailed description of 
the connection between these two definitions and the direct microscopic approach with the purpose 
of serving as a brief referential analysis of this very important method. 

The general expression of the inverse Laplace transform can be written as|19j 



1 l-Pr-ioO 

F{E) = C-^{f{l3)) = / e'^^f{P)d^ , where Pr = real constant , (LI) 

while the part of the expression for the nuclear densities associated with the total nuclear energy 



in (2.11), using the Cauchy theorem and the grand canonical generator function, has the general 
form 

where N is an integer connected with the nuclear energy through E=Ne, defining E as discrete 
with values proportional to a small energy unit, e. Then, one performs the following transformation 
of variables 

y = e-P^ (1.3) 
and assumes y has constant absolute value (constant nuclear temperature) and write 

and, by hypothesis, one has 

/3r.e = logf -j— j- ) = fixed and /3je = — 0„ = — arctgf — ) , (L5) 

Vly|/ \yrj 

then, 

-ed(3 = -eid(3i = idOy . (L6) 

If interval of integration for the argument of y is conveniently taken as 

0^G[-7r,+7r] , (L7) 
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then the integral over the closed path in the y-plane becomes 



(• • ■)dy = I (• • ■)\y\e''y^d9y = - / (• • -Me^'yedP (I.. 



and using (1.4) results 



f-l3r+n/e 

{■■■)dy= / {■■■)e-^^edp. (1.9) 

Now, if e is taken as infinitesimally small the interval of the imaginary integration becomes arbi- 
trarily large, 

Pi = {-Oy/e) G [-vr/e, +7r/e] « (-oo, +oo), (I.IO) 



and, /3j. can be considered as constant (for very small e this would correspond in (1.5) to \y\ very 
close to "1"). Now, if one considers the function in (|I.2[) the expression becomes 



* ^-^dy =:j- _ , ^ ^^^/^ (1.11) 



1 



f{e-^^)e^^ed/3 , (1.12) 
where we have used the relation E=Ne. Then, by definition, when e becomes infinitesimally small 



1 



and the total expression is divided by e, Eq. (1.12) becomes formally equal to C ^(/(e ^^)). For 



(1.12) to be a true inverse Laplace transform one must still assume that E G (0,cxd) so that the term 
proportional to y^ in (2.17) will also be the Laplace transform of the nuclear density for given A 



and E. 

For a given microstate of the grand canonical ensemble, with given finite total energy E 
and total number of excitons n, the above argumentation would imply that becomes very large 
when e becomes infinitesimal. This idea corresponds precisely with the idea of the continuum 
approximation {CAP), and there is no conflict with the direct microscopic description. 

A possible problem in the above deduction could be that /3r= (— l/e)log(|y|), then if e becomes 
infinitesimally small, (3^ would become very large unless log(|y|) becomes infinitesimal too. Now, 
if one interprets as 1/kT, then /3,.— )-oo would mean T— )• 0, which is a physically non necessary 
condition for CAP. 

On the other hand, one may also think that this is irrelevant for our present discussion because 
the variable of integration is /? instead of /3e, therefore, only what affects /3 directly could affect the 
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general argumentation. If log|y| is assumed to be proportional to e with finite and constant f3r this 
would mean that the correspondint points in the y-plane belong to a circle infinitesimally close to 
"1", without essential contradiction. 

The whole reasoning supporting this "interpretation" is therefore coherent with the direct 
microscopic description. 



1.1 The Saddle Point problem 



We saw in equations (2.10) or (2.12) that in formulation of the Shell Model based on the 
Darwin- Fowler statistics the grand canonical generating functional f{x,y) can be written as a sum 
of terms of the type {xy^^), where Vi is an integer. In references [3] and [20] one finds that these 
terms are supposed to decrease in modulus when x or y vary over complex circles around the origin 
in comparison with their values when they belong to the positive real axis and, therefore, they 
should have a maximum in this direction. This assumption corresponds to the idea of a saddle 
point located on the positive real axes of the variables x and y, in the original analysis of Darwin 
and Fowler. l20l 



The formalism of Darwin and Fowler [20] considers the existence of a saddle point for the 



integrand of (2.11) in connection with the steepest descent method. It assumes the existence of a 
point of minimum along the positive real axis and, using a qualitative analysis of the generating 
function of the ensemble, concludes that this point is also a "strong maximum" along the direction 
of the path of integration, taken to be the circle centered at the origin with radius equal to the 
abscissa on the positive real axis where the minimum occurs. 

Reference |20j then concludes that if the integrand has a saddle point with these character- 
istics, it would be sufficient to obtain approximate equations of state in the usual form, for given 
A and E', in the case of systems obeying the Bose-Einstein statistics ("Planck vibrators"). 

Here we briefly analyze the assertion of the existence of this "strong maximum" in the case 



of systems of many fermions, corresponding to the f{x,y) given by (2.10) or (2.12), and show that 
although the saddle point may exist the maximum cannot be characterized as "strong", in the 
sense of having much larger magnitude than the other points along the path of integration. The 
fact that the approximated equation of state obtained with the Darwin Fowler method works very 
well in various applications should then be the object of an independent analysis, which is beyond 
the scope of the present paper. [T8] 

A saddle point of a function is one that is stationary but not a local extremum. For functions 
of two variables, it is a maximum for the variation of one variable and a minimum for the variation 
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of the other. More precisely, it is a point (x^,?/*) G T^"^'" that satisfies [2T] 

L{x^,y)<L{x^,y^)<L{x,y^) , Vx e 7^'^ and Vy E 7^^ (1.13) 

or, equivalently, 

mm^maXyL^x, y) = L{x^, y*) = maXyin.m^L{x, y) , (1.14) 

and the definition for complex variables would correspond to the dimensions n=m=2. 

In a first stage of the formal application of the Darwin-Fowler method the variables x and y 
can be thought as possessing no direct physical meaning and are created only to keep track of the 
counting of the number of particles and the energy of the nuclear levels. This tracking is formally 
performed using the Residue's Theorem in the definition of the "density" associated with y for a 



given A (number of particles, excitons, etc.). Therefore, expression (2.12) for the nuclear density 
would be exact if x and y could be considered as continuous variables, but we have seen that this 
is true only in the CAP limit. pO] 

In addition, the Darwin-Fowler method go one step ahead and attempts to obtain a simpler, 
more practical algebraic expression by performing a qualitative analysis of f{x,y) and using the 
"steepest descent" method, which we will show now can only be considered as approximately 
correct. 

In the following we perform a quick numerical analysis of f{x,y) to show that in at least one 
particular case the hypothesis of the existence of a saddle point with strong maximum along the 
path of integration at the point located on the positive axis is incorrect and, therefore, it cannot 
be considered as having the general validity suggested by Ref.[20j. One can analyze the general 
expression of f{x, y) by considering a simplified version of it, for example, a product of the type 

/ N fix,y) ^ffil + xy'') 

r{x,y) = =11 (1.15 

xy ^-^ xy 

which would correspond to a simplified "ensemble" with maximum number of excitons per config- 
uration equal to "k^ax" and exciton energies ek=ke, with fixed e. The division by "xy" corresponds 
to the description of a system with no excitons, i. e., the fundamental state. We then consider the 
variables x and y as belonging to circumferences around the origin in their respective complex planes 
and variable radii for these circumferences to have the analysis extended to their corresponding 
circles. 

This function can be easily calculated numerically by fixing one variable, for example fixed |x| 
and arg(x), and varying \ y\ with values smaller that one (in accordance to its physical interpretation 
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as y=e~^'^, /3 and e real) and variable arg(y) in a reasonable range, for example between and 27r, 
with steps of 27r/16. 

In figure [s] we show a typical result for A;niax=20 and arg(j;)=0, in which there are indeed 
points of minimum along the real axes of x and y, i. e., for arguments of x and y equal to 0, 27r, 
etc.. In addition, the points of the curve along the positive axis are also the maxima with respect to 
the domain of points situated along the circumferences centered at the origin to which they belong. 
Therefore, these are saddle points in the respective complex planes. On the other hand, one can 
see that these maxima are not "strong" in the sense that one cannot assume that the entire path 
integral along these circles can be reduced to the contribution of a small region in the vicinity of 
these points. Therefore the reasoning of Darwin and Fowler is only approximately valid in this 
case. 
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Figure 8. Simplified generating function, r{x,y), for variable y and fixed x, with arg(x)=0, in the 
case of 20 sp-states. The lines in the box on the right side of the graph are the values of arg(y). 



These conclusions are more clearly shown in figures [9] and 11 below. Figure [9] is the 3 
dimensional correspondent of the previous graph. 

In figure [s] the curves for arg(y) equal zero, 27r, etc. have a minimum for some point between 
zero and 1 and if x is real, this point is also the maximum along the circumference around the 
origin to which it belongs. Therefore, this minimum is a saddle point. 



Notice that the maxima along the real axis in Fig. 8 have not much larger magnitude than 
the other points along the corresponding circumferences, less than 10 times greater for \y\ < 0.7, 
in the given example. The relative magnitude increase for \y\ close to 1, corresponding to nuclear 
temperatures tending to oo, but these points are not saddle points with respect to y. 
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Figures 10 and 11 are the analogous of the previous two figures, now with arg(x)=7r. Notice 
that the maxima along the real axis in comparison with other points of the corresponding circum- 
ferences have completely disappeared and they have become minima instead. In addition, there 
are no points of minimum along the real axis anymore, or in fact in any other direction of the y 
plane. Therefore, for arg(a;)=7r the reasoning of Darwin- Fowler method becomes totally incorrect. 

Therefore, in general, although the maxima on the real axis of x and y may be saddle points 
of r{x,y) they cannot be considered as strong maxima and, if x is considered a complex variable, 
the saddle points can only be defined in certain directions of its complex plane. 




Figure 9. The three dimensional correspondent of Fig. 8 for 20 sp-states with arg(x)=0. Notice 
that the real axis is clearly the loci of local maxima in comparison with the other points along a 
circumference with fixed absolute value of y 
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Absolute value of function JT|^^^^°(1 + x y'^), where both x and y are complex 
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Figure 10. Function r{x,y) in the case of 20 sp-states, with arg(x)=-7r. The hnes in the box on 
the right side of the graph correspond to the values of arg(?/). 
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Figure 11. Three dimensional correspondent of Fig. 10 for 20 sp-states with arg(x)=7r. Notice 
that the points of the real axis are not anymore the loci of local maxima along the circumferences 
with fixed absolute value of y. 



This simple analysis indicates that the entire idea of using the Cauchy theorem in the analysis 
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of the generating function of the grand canonical ensemble, and the consequent connection with 
the Laplace transform, can be misleading. In this case, the use of more direct algebraic approaches 
as the one presented in this paper is naturally more appropriate for a microscopic description of 
the Shell Model, even in the CAP limit. 
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Appendix II - The basic convolution relation 



The expressions for the degneracies in this section result from the combinatorial analysis of 
the corresponding sets of configurations. They are based on the assumption that the maximum 
number of sp-states for particles and holes in the Model Space is arbitrarily defined and the actual 
number of excited particles and holes used to describe the nuclear system corresponds to the 
simple combinatorial idea of placement of "balls" (effectively excited particle or hole) into "boxes" 
(available sp-states for excitation). 

From the operators defined in (|5.4[) it was shown that (a|a^) is given by Eq.(5.33) 



{a\a^) = ^ '^{ph\a'^ai3\p'h'\p'h'\a'^aa\ph) (II.l) 

(12) (1/ A/) (a/3) 

~5^5^e[^^^1^0„^d(p-l + l/a = /3,/i,5,t) , (II.2) 

(1) {S,t) (a/3) 

with S={Ua-{f-a-£i3)) and f=(MQ,-(mQ,-ni/3)). This expression can be rewritten more explicitly using 
the identity 

d{p-a + 13, h, U(^a,^p) -ea + ei3, Mi^^^^p) - + m/3) = 
d{p - a//3, /i, U(^a,^p) - Ea, M(«_-^^) - m„) = 

d{p - a, /i, Ua - Ca, Ma - hIq,) - d{p - a = h, Uap - ea^M^p - m„) , (II.3) 

where (C^(Q,^/3),Af(^ -^^)) are the energies and angular momenta of the configurations in which (a) 
is present and (/3) is not and {Uap^Map) are the parameters of the configurations in which both (a) 
and (/?) are present, etc. If one considers the set of all configurations in which (a) is present they 
have degeneracy, d{p=a,h,Ua,Ma), equal to the first term of the last expression of ( |II.3 ) 



d{p = a, h, Ua, Ma) = d{p — a, h, Ua — Ca, Ma — IHo) (II-4) 

and it includes also configurations in which (/3) is possibly present. The degeneracy of the config- 
urations in which both (a) and {(3) are simultaneously present is 

d{p = a(3, h, Uaf3, Mali) = d{p - a = (3, h, Uafi - Eq,, Map - m^) (II.5) 
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therefore, the subtraction of (II. 5) from (II. 4) produces the degeneracies of the configurations in 



which (/3) is not present, which is (II. 3). Note that one can also write (II.3) as 



d{p = a//3, h, Ua^i3, Ma^p) = d{p = a, /i, Ua, Ma) - d{p = a/3, h, Uap, Ma 
then ( |II.2[ ) can be rewritten as 



(11.6) 



e^^^'^^'^Oa(3[d{p - a,h,Ua - ea,Ma -niQ,) - d{p - a = /3,h,Uai3 - ea,Ma(3 -m„)] , (II.7) 

\{S,t)) 



and, if Oa(5 is given for ah (II.7) can be determined from the combinatorial results for 



d{p-l,h,U-ei,M-m.i) and d{p-2,h,U-ei-e2,M-mi-m2) 



Similarly, the corresponding term for holes, (5.39), can be expressed as 



(e|e+)= T.(pf'H^-\P'''lph'\bt.bi3\ph) 

{12)iUM)(a(S) 

- ^ e[^*^1^0„^d(p,/i-l + l//3 = a,5,t) , 

( ) (a/3) 
\(S,t)) 

with S={Uis-ei3+ea) and t=(M^-m^+mQ,), which can be rewritten as 



(11.8) 
(11.9) 



el^^^^YOMdiP, h-P,Up- ep,Mp - m^) - d{jp, h-P = a,Ua(s- e/3, - m^)] , (11.10) 

\{S,t)) 

The above expressions can be rewritten as convolutions with the various sp-state densities, 
for example, if one consider the following "interpretation" of the sum over one sp-state. 



^/(ea,ma) = J c^e^^5(e - e„)/(e, 

(a) (inQ,)(eci) 



m 



(11.11) 



and 



where usually the sp-state "a" is degenerate with respect to but not with respect to the pair 
(eQ,,niQ,). Then, 

^/(ea,ma) = j (ie^/(e,ma)cj(l,0,e,ma) , (11.13) 

(a) (m„) 

with the following straightforward inclusion of a second sp-state, 

^/(ea,e/3,ma,m/3) = j deide2 ^ /(ei, £2, mc,,in^)u;(l, 0, ei, ma)a;(l, 0, £2, m^) , (11-14) 



(a/3) 
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where all sp-energies are positive and the sp-densities are functionals that restrict the contributions 
of the degeneracies to points where they are well defined. 



Then, one can rewrite the sum over {a, (3) in the first term of (II. 7) as 

'^Oafsd{p - a, h, Ua-€a, Mq 



' a ) ) 



(11.15) 



{al3) 



= jdeide2 ^ Oai30j{p - i,h,U - €i,M - ma)uj{l,0,ei,ma)uj{l,0, €2,^13) , (11.16) 

(mc«m^) 

j ^^iX^f j "^^2^^(1,0, €2, ni;3)Oc«/3 |'^(p- l,h,U - ei,M -ma)u){l,0,ei,ma) , (11.17) 

° (m„) \ (m^) / 

r-U 

= deiJ^OaCoip - hh,U -ei,M -ma)uj{l,0,ei,ma) , (11.18) 

Jo , ^ 



where 



(ma) 



a;(l, 0, t2-,'^^)0ap I , 



(11.19) 



and the sum over (a,/3) in the second term of (II. 7) as 

Jdeide2 ^ Oa/3w(p - 2, /i, ^7 - ei - e2, M - nia - m/3)a;(l, 0, ei, ma)tj(l, 0, e2, m^j) (11.20) 



U 

JO 



del (ie2 Oa/3a;(p — 2, /i, [7 — e, M — ni)u;(l, 0, e — £2, m — ni^)a;(l, 0, 62, m^) (11-21) 



= / deVOe(2,0,e,m)w(p-2,/i,C/-e,M-m) , 

(m) 

where e=(ei+e2) and ni=(mQ +111/3) and 

O,(2,0,e,m) = / (ie2 ^^©^/^^^(l, 0, e - €2, m - m/3)w(l, 0, e2,m;3) , 
JO 



(11.22) 



(11.23) 



(ma) 



Therefore, both terms can be expressed as convolutions if adequate definitions are given for 
the nuclear degeneracies. 

In |Sec.4 a different approach is used for the nuclear density, without Dirac deltas, which is 
equivalent to the above one. For each (U,M), the density was defined by the following approximated 
expression. 
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U, M) w d{A, U, M)/6U = constant, 
which imphes that the integral of uj{A,U ,M) over [Uprev,U] is equal to the degeneracy at U, 

d{A,U,M) M 

d{A,U,M)= (1) ^uj{A,U,M)6U = uj{A,U,M)dU. (11.24) 



then, 



Ema. d{A,U,M) E^^^ 

E (•••)^ / u:{A,U,M){...)dU. 

{U=E^,„) (i=l) 



(11.25) 



U=E„ 



In the case of sp-states a similar definition can be used in which the degeneracy in (11.23), at 
each point (ec,nia), is usually equal one, then 



1 



w(l, 0, e, ma)de. 



(11.26) 



where uj{l,0,e,ma)=l/5e, and 6e={ea-ea,prev) and 



(a) (eci=e7nin){ma) 



y^a;(l,0, e,nia)(...)de 

(ma) 



(11.27) 



Then, replacing (...) by ^^Oad{p-a,h,Ua-£a,Ma-may\ with Oa given by (11.19), and using (11.24) 
yields 



Y^adip - a, h, Ua- ea, M, 

(a) 



a-ma)= Y ^Oad(P - a, ^, f^a - Afa - ma) (11.28) 



^^a;(l, 0, e, ma)Oa / uj{p — a, h,U — e, Ma — ma)dUd€ , 

-min ^ Ua^prev 



(11.29) 



and if one assumes that the approximate densities vanish outside the energy interval used in their 



definition, one can rewrite (|II.29|) as 

'U. 



/ dU^Oa a;(l,0,e ^ma)uj{p — 1, /i, J7 — 6, Ma — nia)de 



(11.30) 



and if one assumes that the convolution integral is approximately constant in the interval {Ua,prev,Ua\ 
results 



YiOadaip - a, h, Ua - M„ - m^) 
(a) 



(11.31) 
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l-u 

K.5UY]Oal u{p-l,h,U - e,M -m)Lo{l,0,€,m)d€ (11.32) 

(ma) 

where 5U ={Ua-Ua,prev) 1 then the sum of "^(q,)©^*^^" is approximatctly proportional to the con- 
volution of nuclear densities. Now, under CAP one may assume that 5U is approximately constant 
and small, 

5U 

~ ^nuc — infinitesimal nuclear energy spacing (11.33) 
then, taking 5U ^ dU results 

J^Oada^^^ i^^Oa f uj{p - 1, h,U - e, M - ma)cvil,0, e,mo,)de\ dU (11.34) 

{«) \{m„) J 



which can be interpreted as the elementary variation of the inverse Laplace transform of the summed 
transition strengths. 
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Appendix III - Analysis of the Model Space description 



In this section the expression of the one body operator 

0= Oapaiap (III.l) 

/a e v\ 

is analyzed and its Model Space interpretation, based on the combinatorial analysis of the occupa- 
tion of sp-states, is compared with the physical interpretation of the operators. 

To obtain an expression involving also the hole operators so that transitions in which "new 
particle" and "new hole" states can be excited, as opposed to transitions corresponding to the 
simple scattering of the already excited "particle" states, one introduces the operators associated 
with the "/i" states and replaces 

by (a^ + ftj^) and by {al + 6„) (111.2) 



then, 



O = ^ Oai3aiai3 + ^ Oai3aib^i3 + ^ O^phaap + ^ Oa/sbabl ■ (III.3) 



(aep\ (aep\ (a£h\ { a £ h\ 

\l3epJ \i3eh) \fiep) \iseh) 



When compared with (III.l) the last expression needs some clarification because, on the one 



hand, the statement (aGp,/3Gp) in (III.l) means that a "p" was destroyed in (/3) and another "p" 
was re-created in (q), therefore both and (q) are sp-states of the "p" type in the Model Space. 
On the other hand, if (a) is not equal to (/3), does the possible creation of "p" in (a) necessarily 
means that (a) is "empty" in the sense that that it is a "/i" sp-state? 

To answer this question one may think initially within a combinatorial point of view and 
consider the "p"-type states defined by the Model Space as "empty boxes" in which the various 
"p" excitons will be placed when the Oq operator acts and, similarly, the action of the operator 
will decrease the number "p" excitons but will not alter the number of "p"-type states originally 
defined in the Model Space, which is also the maximum number of "p" excitons in the model. 

The interpretation of ba and ba in the Model Space is similar to the "p" operators with the 
caveat that one is dealing now with the placement or withdrawal of "holes" into and out of "empty 
boxes" , which may not be a physically intuitive idea. 

In a physical point of view a "p" sp-state can only be created on states occupied by "/i" 
sp-state and vice versa. Then, when ba operates it can act exclusively on states occupied by "p" 
although in combinatorial terms the action of ba means only that a new "/i" state (a) has been 
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"placed" into one of the available "empty boxes" for holes of the Model Space. Then, which one is 
correct? To say that the domain of sp-states on which ba acts is "a Sp" or "a E /i"? 

In the Model Space description "particles" and "holes" are considered independent fields, 
then the creation of "/i" can happen wherever there are "empty boxes" and the consideration of its 
"previous occupation by a particle" becomes meaningless. Then, one can use the fermion relation 

bibfs + bf3bi = 6a,p , (III.4) 



where both (a) and (/3) are "/i-boxes" of the Model Space ("empty" or "filled"), with no direct 
relation with the "p" sp-states. The case a=f3 corresponds either to the action of "filling a box 
that has just be emptied" or "emptying a box that has just been filled" and the case a^P gives 
the "zero operator" meaning that the both processes of creating a "/i" (a) followed or preceded by 
the destruction of a different "/i" (/3) have zero probability. 

On the other hand, if one thinks more physically that ba is acting on states occupied by "p" 
and ba corresponds in fact to the creation of a "p" on states previously occupied by holes, then the 



first term of (III. 4) is acting on the "/i" sp-states that have already been excited, while the second 
term is acting on the "p" sp-states that also have already been excited and neither of them are 
acting over the entire set of "/i-boxes" of the Model Space. 

In combinatorial terms the action of on \ph) selects configurations (p) that contain (a) 
while ba acting on \ph) selects the configurations (h) that do not contain (a) . Then, (aa+ba) 
acting on \ph) select the configurations (p/i) in which (a) is present in (p) or not present in (/i) 
and the number of such configurations is 



e::)(G)-G!0)^C:O(^)(::O^e::)Cr)' 

where the maximum number of "p" excitons in the Model Space is being called g and the maximum 
number of "/i" excitons is being called b and these letters will be used exclusively to designate these 
numbers in this paper. These maxima are related with the physical limits for the possible energies 
of excitons beyond which the excited "p" is supposed to be emitted and the nuclear system will be 
essentially altered in the sense that the Model Space must be modified too. 

Because the sp-operators Oq, and ba act on different fields and the fact that the action of Oq, 
decreases the excitation energy, U, while ba increases U, one concludes that the sp-states designated 
by (a) in both operators are in fact independent, therefore the use of the same index (q) to designate 
both sp-states can be misleading. 
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The dynamics of the whole many-body system is defined independently of the Model Space 
itself. For example, the idea that "p"s and "/i"s can only be created or destroyed in p/i-pairs 
is described in the model strictly by additional constraints to correctly describe the conservation 
laws and the change of the number of excitons, with the corresponding change of U and angular 
momenta, by the "correct definition" of the residual interaction. 

Then, the (a) on Oq acts on the "p" sp-states which have already been excited in the available 
nuclear configurations {ph) while the (a) of ba acts on the "/i" sp-states of the Model Space that 
have not yet been occupied by a "/i". Similarly, in the case of {a)a+ba), Oa acts on the "p" sp-states 
that have not been excited yet while ba destroys the already existing "/i" . 

Now, when one puts both operators together 

{ai + ba){aa + bi) , (III.6) 

the interpretation must follow logically from the above one. In particular, the term corresponds 
to the creation of a "/i" on an unoccupied state (a) of the Model Space followed by the destruction 
of the same sp-state. Therefore, it counts the number of unoccupied ( "available" ) sp-states of the 
type "/i" while h\j)a counts the number of occupied sp-states of the type "/i". The sum of these 
operators, (6L&a+&a&a), counts the total number of sp-states of the type "/i" in the Model Space 



and, using the notation of (III. 5), this number is b. 

The fact that the excitation of "p"s and "/i"s is physically not indenpendent of each other 
complicates the description a little because b is usually greater than g (because one cannot define 
more "p" sp-states than the actual number of nucleons in the nuclear system, while the number of 
"/i" sp-states has no constraints except for the maximum nuclear excitation energy) and, on the 
other hand, for a given Model Space, the maximum excitation of "/i" sp-states cannot in fact go 
beyond g. Then, 



(ph\b%\pli) + Y,<P''\hl'i\p'>) = J2 (pl'\Sa.f\ph) = b , (III.7) 



(a/3 G f illed-/i) (a/Se-h) (a/3Gall-h) 



and the normal ordered operator O in Eq. (III.l) can be rewritten as 



( ae-p \ (aefilled-p) (aG-p\ /a S f illed-hN 



- ^ Oal3b^pba - 



In reference to (III. 8), the last term is 
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Oaa= ^ {a\Oal3Sa,l3\(3) , 
(«eall-/i) (a,/3 e all-/i) 



and if one considers, for example, the number of particles operator (III. 8) gives 
N = 



then 



(a e f illed-p) (a S f illed-h) (aSall-h) 



(AT) = (p/i|iV|p/i) =p-h + b, 



(III.9) 



(III.IO) 



which agrees with the idea expressed in (III. 5) and is equal to the total number of already excited 



"p" sp-states ("filled" p states) plus the total number of "/i" sp-states that have not yet been 
excited ("empty" h states). This description equates filled "p" sp-states with empty "/i" sp-states, 
in which "/i" sp-states can be created. This would be physically correct if "/i" sp-states could only 
be created by the destruction of a "p" in the same sp-state, but this is not necessarily true in the 
Model Space description. 

In other words, the total number of actually excited "p"s, which is the expected value of 
Qaaa, has been increased by the total number of "potential holes", which coincides in the Model 
Space with the total number of "empty-/i" sp-states, but this reasoning does not take into account 
the constraint due to pair creation and annihilation. If this condition is taken into account, the 
actual number of potential holes ( "empty-/i" states) cannot be greater than "g-p" instead of "&-/i" , 



corresponding to the interpretation of holes as empty particle states, then (III.IO) becomes {N)=g, 
which is the physically correct result. 

In conclusion, the apparent simplicity of the basic formulation of the Model Space that leads 



from (III.l) to (III.3), may produce meaningless results if it is not epistemologically enhanced by 



the correct physical interpretation of each step. This happens due to the inherent incompleteness 
of the Model Space description. 
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Appendix IV - A problem with the Laplace transform approach 



In the following discussion we drop the dependence of the degeneracies and densities on the 
total angular momentum, M, to have a clearer and easier to follow presentation, but it could be 
included in the argument in a straightforward way. The following analysis was inspired by the 



results for the evaluation of the moments of the Hamiltonian, in Sec. 6, in which a given sp-state is 



destroyed and created in the same transition, in particular by Eqs. (6.17), (6.18) and (6.19) 



In general one may write the generating function of the grand canonical ensemble as 



(1,1/ 



(p,h=0) (ei<--<e(p+,i)) 



(IV.l) 



where U=ei+- ■ •+e(^p^h)i and for a given "particle" sp-state (a) one may write 



fa{x,y) 



1 + X 



pa 



oo 

E- 

{p,h=0) 



(IV.2) 



where U is the energy of p-particles and /i-holes with all sp-states considered, while U^^a is the 
energy of p-particles and /i-holes in which all p-particles are different from a. Then, one can 



rewrite (IV.l) as 



f{x, y) = (1 + Xpa)fa{x, y) = fa{x, y) + Xpafa{x, y) 



(IV.3) 



and the last term corresponds to the terms summed on the RHS of (6.17). 



Notice that for given the energy grids of U-^a and Ua are "complementary" with respect 
to the grid of U in the sense that if g and h are the maximum values of p and h in the Model Space, 
then one can show by combinatorial analysis of the corresponding sets of configurations that the 
number of configurations associated with the domains of [/, Ua and C/^q,) which will be called 
NOC{U), NOC{Ua) and NOC{U^a), respectively, are given by 



NOC{U) 



, NOC{Ua 



g-1 



and NOC{U^c 



g-1 

P 



\Pj\hJ \p—lj\hj \ p J \h 

where the conventional notation for binomial coefficients has been used. Therefore, 

NOC{U) = NOC{Ua) + NOC{U^a) , 



, (IV.4) 



(IV.5) 



and NOC{U) is g times greater than NOC{U^a)^ corresponding to very different grids usually. 



From (IV.2) one derives an approximate connection with the Laplace transform 



^.e-P^-fa{x,y)= Y ^"^V d(p/a,/i,C/^«)e-^(^^"+^'^) 

(p,h=0) {e^^e^) 



(IV.6) 
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CAP 



x"" y 



ip,h=0) 



poo 

Jo 



'-a, h,s — €a)e ^'^ds 



(p,h=0) 



(IV.7) 
(IV.8) 
(IV.9) 



(p,h=0) 

where "ti(s-eQ,)" is the Heaviside step function and "aj(p^a,/i,i?)" is the continuous nuclear density 
corresponding to the discrete degeneracy "d(p/a,/i,C/^Q)" (and "ii" is the continuous counterpart 



(IV. 10) 



a, as defined in |Sec.4| ). On the other hand, from ( |6.18[ ) one may write 

Xpafa{x,y) = xe-^'-Ux,y) ''^''^^e^^'^e-^'-d{p-l^a,h,S,t) 

(1) {St) 



^c{e-^^-u;{p-l^a,h,S,t)] , 
(1) 



(IV.ll) 



where S={Ua-€a)- This expression would again lead to a Laplace transform involving the Heaviside 

(IV.12) 



step function, but one could also use (6.19) to obtain 

CAP 



Xpafaix,y) = xe ^^"/^(x,?/) ^C{u;{p = a,h,Ua,Ma)} 

(1) 



which, due to the identity (6.16), is incompatible with (IV.ll) unless ea=0. 



Equation (IV.12) is the correct one in this case because it is in agreement with (IV. 5) while 



(IV.ll) is not, as it will be demonstrated next. 



Initially notice that in (IV.12) the nuclear density can also be written as u>{p-a+a,h,Ua,Ma) 
because the selection of configurations is determined by the destruction of (a) and not altered by 
the following re-creation of (a). 

On the other hand, the degeneracy "d(p7^a,/i,[/^ci)" is not the same as ^^d{p-l^a,h,Ua- 
Ca)" in (6.18) because the energy grids of the arguments of both functions are different (in fact, 
complementary in the sense of ( |IV.4 ) and (IV. 5)), but it is easy to show by combinatorial analysis 
that "d(p7^a,/i,[/^a)" is equal to "(i(p-|-l-a,/i,[/^a)" i corresponding to the grid of energies obtained 
by the creation of a new particle in a given configuration, increasing "p" to "p+l", followed by the 
"selection" of configurations containing (a) by the destruction of (a). 

From (IV.9) one has for given (a) and {p,h) the following inverse Laplace transform of (IV.4) 

C~^[xe-^'-Ux,y)} = C-'{e-f''''CMp^a,h,U^a}} =u{U^aMp^a,h,U^^^ , (IV.13) 
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and using £ ^{fa{x,y)}=uj{py^a,h,U:^a), ^ ^{f{x,y)}=uj{p,h,U) and (IV. 5) one would conclude 
that 



uj{p, h, U) = uj{pi^a, h, U^a) + u{U^a)^{p^a, h, U^a) , 



(IV.14) 



which is obviously a false relation because (q) is present in the domain of the configurations 



described by the LHS but is not present on the RHS. On the other hand, according to (IV. 5) and 



using (IV. 12), the correct relation would be 



uj{p, h, U) « Lo{p - 1/a, /i, Ua - to) + ^{pi^a-, h, U^a) 



(IV. 15) 



or 



uj{p, h, U) w uj{p = a, /i, Ua) + (^{pT^a, h, U^a) , 



(IV.16) 



which, in this case, is an obviously true relation. 

This result means that the exponential term involving the sp-energy, in the term Xpafa{x,y) 



of (IV.7), should be "absorbed" into the integral defining the Laplace transform instead of be- 



ing considered as a multiplicative term that produces an inverse Laplace transform in which the 



Heaviside step function is present. This result corresponds to the identity between Eqs.(6.18) and 



(6.19). 



Therefore, the direct use of the Laplace transform in the analysis of (a) in (6.12) may lead 
to wrong conclusions. 

Sometimes the creation of sp-states brings problems of interpretation under the Laplace 
transform approach, which can be easily overcome within the direct microscopic formalism. 
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Appendix V - Transition strengths expressed as convolutions 



This Appendix is directly related to Sec.6 and presents some details of the algebraic deduc- 
tions of a few transition strengths (TS) of the PE Hamiltonian, in particular, of how they can be 
expressed as convolutions with nuclear densities of less excited states. Each TS is preceded by a 



title identical to the corresponding subsection of Sec.6 



Terms that increases the number of p and h by 2. 



The expression of (a|a ) in (|6.61 ) is 



(a|a+) = ^ ^ {ph\a^aphjh1:^\p'h')j)h'\h.ybsai3aa\ph) 

(-,)M<57) 

= EE^'""*'^ E \yc.M&\^d{p - 2/a/3, h - 2/^7, 5, t) , 

(1) (St) (a/357) 



(V.l) 
(V.2) 



which are analogous of the first moment expressions. 

This expression can be rewritten, using (??), in terms of the convolution over the energies 
of the selected sp-states by rewriting the sum over {aj3"f5) as. 



Yl \VaMS?d{p-2^aP,h-2^5^,S,t) 
(0/3-57) 

^ \Vap-i6\'^d{p -a- /3,h-6-J, U(^a,p,-/,5) - e, ^(a,/3,7,5) 
(q/37<5) 

J deide2desde4 \Vai3-ys\'^d{p — 2,h — 2,U — e, M — 

(mcm^m-ym^) 



(V.3) 
(V.4) 



xw(l,0, ei,ma)w(l,0, e2,ni^)tj(0, 1, 63, m^)w(0, 1, €4,1115) 



(V.5) 



= E / deV{2, 2, e, m) u}{p - 2,h - 2,U - e, M - m) 

where e=eQ+e^+e-^+e5, m=nio+ni^+m^+ni5 and 

y(2,2,e,m)= ^€4 | V^^^5p 0, e — x, m — nix) 



(mx ,my ) 



(V.6) 



X ui{l,0,x- 2/,mx - my)u;(0, l,y - e4,niy -ms)uj{0, 1,64,1115) 



(V.7) 
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and y=e3+€4, x=e2+y, mx=ni^+my and niy=m^+mo 



Equation (V.6) cannot be further simplified unless explicit assumptions are made about the 
dependence of |V"a/3^5p on the energy and angular momentum of the sp-states. 

Terms that increases the number of p and h by one, 



From (6.76) one has 



(6|6t) = ^ ^ \Vap^s\'^{ph\a'^a'l^bja^\ph"lp'h'\a'^bsai3aa\ph) 

(-,)(«/3^7) 

= EE^^''''" E \Vap^5\''d{p-a-P + j,h-5,S,t) . 

(1) (5,t) {al35-y) 



and one can rewrite the sum over {aP^6) as (see (III. 3)) 



\Vai3^s\'^ {d{p — a — /3,h — 6, S,t) — d{p — a — /3 — ^, h — 6, S,t)) . 
The first term is (see (??)) 

\Vai3^s\'^d{p - a - f3,h - 5,S,t) , 

J d€ide2de3de4 \Vai3'y5\'^d{p — 2,h — 1,U — e, M — m) , 

xa;(l, 0, ei,mQ,)a;(l, 0, e2,m^)a;(0, 1, e3,m^)a;(0, 1, 64,1115) 
= / d€V^{2,l,x,:a^)ujip-2,h-l,U -x,M -m^) , 



(V.8) 



(V.9) 



(V.IO) 



(V.ll) 



(V.12) 



where 2/=(e2+e4), x={ei+y), nix=(mQ,+my) and my=(ni^+m5), 

14(2, l,x,mx) = y dy d€4V(^ai35) w(l,0,x - y,mx - my) w(0, l,y - £4, my - m5)a;(0, 1, £4, m5)(V.13) 

(my, ma) 

and 

y(ai35) = j (ie3^|Va/3^5pw(l,0,e3,m^) , (V.14) 



and the second term of (V.IO) is 
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\Vai3-ys\'^d{p - a - f3 - -f,h- 6,U - e,M -m) , 

^ J d€ide2de3de4 \Vai3'y5\'^d{p — 3,h — 1,U — e, M — m) , 

xa;(l,0, ei,ma)a;(l, 0,62, ni^)tj(0, 1,63, ni^)w(0,l, €4,1115) , (V.15) 

= Yj deV{2, 2, e, m) d{p - 2, h - 2,U - e, M - m) (V.16) 
W 

where e=(ei+e2+e3+e4), ni=(mQ,+m^ +111^+1115), ^=(£3+64), x={e2+y), nix=(m^+my) and niy=(m^+niQ) 
and 

y(2,2,e,m)= J dx dy de4\Vai3'ys\'^ oj {1,0, e — x,m — m^^) 

(nix ,niy ,ms ) 

X a;(l,0,x - y,mx - my)a;(0, l,y - e4,niy -1115) a;(0, 1,64,1115) (V.17) 

Terms that increase {p,h) by one: off-shell processes. 

In the cases in which Hermitean conjugates of a given sp-operator, apa and Opa, act on the 
same nuclear configuration the corresponding sp-energies are not present in the energy conservation 
constraint for the total transition. Then, the processes of creation and destruction of these sp-states 
are off-shell with respect to the transition process. 

In the expression of (6|6^) in ( |6.65[ ) with 7 replaced by a in (6) and (6+) one has 

(6^^a|6+_^„) = ^ {ph\a^a^bjaa\ph'\ph'\a^bsai3aa\ph) 



12 



)(a/35) 



^^^[c/M] ^ \v^^s\^d{p -a-P + a,h-6,S,t) , (V.18) 

(1) {S,t) {apS) 



as it was shown in (6.96), where 



S = {Ua-ep- 65) and t = (M„ - - m^) , (V.19) 
d{p - a - /3 + a,h - 6,S,t) = d{p - a - j3,h - 5, 11(^^/35) - x-, ^(aps) - i^ix) , (V.20) 
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and 

X = ea + ep + and nix = ma + m^g + . (V-21) 

The indices of U and M correspond to the sp-states that are necessarily present in the nuclear 
configurations associated with these parameters, for the transition to be possible. Although the 
configurations after the transition have 2 excitons less than the initial one, the degeneracy describing 
the transition strength corresponds to configurations with 3 excitons less because (a) was selected 
by an additional off-shell process. 

Then, 

(6^^„|6+_^„) « jdxY^V:ri2,l,x,my,)u;{p-2,h-l,U -x,M -iTiy,) (V.22) 



where 



Vx{2, l,x,mx) = J dydes ^ | VQ,^^5pa;(l, 0, x - y,mx - niy) 



(myms) 

xuj{l,0,y - e3,my - m5)a;(0, 1, €3,1115) , (V.23) 



which is expected to be valid in the CAP limit. 



Similarly, another off-shell process can be obtained from (6.85) with a replaced by a different 
index 7 in {b+), i.e., j^a, then 

= ^ ^ {ph\aia^pblaa\ph'Xph'\a\bsapa^\ph). (V.24) 
^EE^"''''^ d{p-a-^-l3 + a + ^,h-5,S,t) , (V.25) 

(1) {S,t) {al3S-y) 



where S and t are still given by (V.19) and now both a and 7 transitions are described as off-shell. 
The degeneracy is 

d{p — a — ^ — 13 + a + ^, h — 6, S,t) = d{p — a — f3 — 'jjh — 5,U — €,M — m) , 



= d{{p = a-i) - P,h- 5, U(^af) - e/3 - e^, M(c,^) -m/s - mg) , 



(V.26) 
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and 

X = ea + + and nix = nia + + . (V.27) 

Although the configurations after the transition have 2 excitons less than the initial one, the 
degeneracy describing the transition strength corresponds to configurations with 4 excitons less 
because (a) and (7) were selected by an additional off-shell processes. 

Then, 

{b^^a\hU^)^^^- y"a!eJ^y,(3,l,e,mV(p-3,/i-l,C/-e,M-m) (V.28) 



where 



(m) 



V;(3, l,e,m) = j dxdyde^ ^ | Va/3^5|^a;(l, 0, e - x, m - nix) 

/raxmy \ 

xa;(l, 0, X - y, nix - my) a;(l, 0, y - £4, my - m^) a;(0, 1, £4, m^) , (V.29) 
which is valid in the CAP limit. 

Other terms with on-shell transitions. 



The term corresponding to (c|c+) in (6.10) can be written analogously 



(c|c+)= J] Y,\V^p^8\Hph\aib\h\hp\p'h'lp'h'\h^^h^ha^\ph) (V.30) 
~5Z^e[^^^l ^ |y„^^^|2d(p_a,/i-<5-7 + /3,5,t) , 

(1) (5,t) (a/3<57) 

^^et^^-^l ^ (|V^^^5|2d(p_Q-,/i-(5-7,f/(„5^) -x,M(c,5^) -nix) 

(1) (5,t) (a/3<57) 

-|Va/3^5pd(p-a,/i-(5-7-/3,[/(„;3^5) -e,M(^ap-i5) -m)) , (V.31) 



where 



5 = ({/ — ea — — £5 + e^) and t = (M — — — + m^) , (V.32) 
e = (ea + 6/3 + + e^) and m = (niQ, + + + m^) , (V.33) 

a; = (cq + + e^) and nix = (mo- + m^ + m^) , (V.34) 

and the indices of U and M indicate the sp-states which are supposed to be present in the nuclear 
configurations associated with these parameters. 
The first term is 
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X] l^a/37'5pc?(P - a,h - 6 -j,S,t) , 

(0/3-57) 

J deid€2desde4 ^ |Fa;375|^w(p - 1, - 2, J7 - x, M - nix) , 

xa;(l,0,ei,ma)a;(0,l,e2,nLy)w(0, 1,63, m5)a;(l, 0,64, ni;3) , (V.35) 

= dxl4(l,2,a;,mx)a;(p- l,/i-2,C/-a;,M-mx) (V.36) 

(mx) 



where 



T4(l,2,a;,mx)=^ / dydes V'(c,^5) a;(l, 0, x - y, nix - my) a;(0, 1, y - 63, my - 1115) a;(0, 1, 63,1115), (V.37) 

(mym^J 

y = {e^ + €s) and my = (m^ + m^) , (V.38) 



and 

and the second term is 



(m) 

where 



X] \'^a^'ys\'^d{p - a,h- P - § -j,S,t) , 
J deide2d€2d€4^ |Kk;37<5|^'^(p — 1, /i — 3, i7 — e, M — m) , 

xa;(l,0,ei,mQ,)w(0, l,e4,m^)a;(0, l,e2,m-j,)a;(0, l,€3,m5) , (V.39) 
J2 JdeVeil, 3,e,m)uj{p-l,h-3,U-e,M-m) (V.40) 
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^(1,3, e,m) = X] / dxdyde^ \Vap^6\^uj{l,Q,e- x^va-my,) 

X 0,x-y,myr- my) uj{0, l,y- £3, my - mj) w(0, 1, €3, m^) 

with (x,mx) redefined as 

x = {ej + es + e^) and mx = (m^ + m^ + m^) , (V.41) 

The terms corresponding to {c\c^) with fi^S are 

{c\c^)p^5 = Yl E (Ph\a+b+b+bp\p'h'lp'h'\b+b^bsa^\ph) 

_ ^ et^^^l ^ - a, /i - (5 - 7 - ^ + /3 + (5, i7 - - e^, M - ma - m^) , 

((l),5,t) (a/357) 
(a/357) 

-\Va0i5?d{p -a,h-6--/- P, Ua^jS - e, Map^s - m)) , (V.42) 



J2 e'^^l y" de^FeCl, 3, e, m)a;(p - 1, /i - 3, C/ - e, M - m) (V.43) 



where 



and 



((l),5,t) (m) 

V;(l,3,e,m) = / dxdyde^ ^ |VQ/3-^5pa;(l, 0, e - x,m - mx) 

xuj{0, 1, a; - y, mx - my) a;(0, l,y - £4, my - m^) a;(0, 1, 64, m^) , (V.44) 

e = ea + + + es and m = m^ + m^ + m^ + m^ , (V.45) 

X = ejs + e-y + es and mx = m^ + m^ + m^ (V.46) 

y = + es and my = m^ + m^ . (V.47) 
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